LOCAL BIFURCATION-BRANCHING ANALYSIS 

OF GLOBAL AND "BLOW-UP" PATTERNS 
FOR A FOURTH-ORDER THIN FILM EQUATION 

P. ALVAREZ-CAUDEVILLA AND VA. GALAKTIONOV 

Abstract. Countable families of global-in-time and blow-up similarity sign-changing 
patterns of the Cauchy problem for the fourth-order thin film equation (TFE-4) 

u t = -V • (|u| n VAu) in l^xli, where n > 0, 

are studied. The similarity solutions are of standard "forward" and "backward" forms 

u±(x,t) = (±t)- a f(y), y = x/{±tf, /3 = ±=p, ±t > 0, where / solve 

(0.1) B±(a,/) = -V-(|/|«VA/)±/3yV/±a/ = in R N , 

and a £ R is a parameter (a "nonlinear eigenvalue"). The sign "+", i.e., t > 0, corre- 
sponds to global asymptotics as t — > +oo, while "— " (t < 0) yields blow-up limits f->0" 
describing possible "micro-scale" (multiple zero) structures of solutions of the PDE. 

To get a countable set of nonlinear pairs {/ 7 , a 7 }, a bifurcation-branching analysis is 
performed by using a homotopy path n — > + in (10. where Bq (a, /) become associated 
with a pair {B, B*} of linear non-self- adjoint operators 

B = -A 2 + \y V+f I and B* = -A 2 - \ y • V (so (B)* 2 = B*), 

which are known to possess a discrete real spectrum, er(B) = u(B*) = {A 7 = ~"^}|7|>o 
(7 is a multiindex in R^). These operators occur after corresponding global and blow-up 
scaling of the classic bi-harmonic equation ut — — A 2 u. This allows us to trace out the 
origin of a countable family of n-branches of nonlinear eigenfunctions by using simple or 
semisimple eigenvalues of the linear operators {B,B*} leading to important properties 
of oscillatory sign-changing nonlinear patterns of the TFE, at least, for small n > 0. 



1. Introduction: TFE-4 and "adjoint" nonlinear eigenvalue problems 

1.1. The model, two classes of similarity solutions, and main problems. In this 
paper, we study global asymptotic behaviour (as t — > +00) and finite-time blow-up be- 
haviour (as t — > T~ < +00) of solutions of the fourth-order semilinear thin film equation 
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(TFE-4) 

(1.1) u t = -V ■ (\u\ n VAu) in R N xl + , n > 0, 

where V = grad x and A = V • V stands for the Laplace operator in Mr. 

Before describing the main application of such a nonlinear higher-order PDE model, 
which has become widely known in the last decades, we state the main goal of the paper. 
We study similarity solutions of (11.11) of two "forward" and Sturm's "backward" types: 

(i) global similarity patterns for t 1, and 

(ii) blow-up similarity ones with the finite-time behaviour as t — > T~ < oo. 

Both classes of such particular solutions of the TFE-4 (11.11) can be written in the joint 
form as follows (here, the blow-up time T = for solutions in (ii)): 



I— an 



(1.2) u±(x,t) = (±t)~ a f(y), y = x/{±tf for ±t>0, where (3- , . 
and similarity profiles f(y) satisfy the following nonlinear eigenvalue problems, resp., 



1.3) (NEP)± : B±(a, /) = -V • (|/|"VA/) ± 0y ■ V/ ± af = 



m 



Here, a G M is a parameter, which stands in both cases for admitted real nonlinear 
eigenvalues. Thus, the sign "+", i.e., t > 0, corresponds to global asymptotic as t — > +oo, 
while "— " (t < 0) yields blow-up limits t — > T = _ describing a "micro-scale" structure of 
the PDE. In fact, the blow-up patterns are assumed to describe the structures of "multiple 
zeros" of solutions of the TFE-4. This idea goes back to Sturm's analysis of solutions of 
the ID heat equation performed in 1836 [37]; see [24"l Ch. 1] for the whole history and 
applications of these fundamental Sturm's ideas and two zero set Theorems. 
Being equipped with proper "boundary conditions at infinity" , namely, 

(1.4) for global case, B^(a, /) : / G Co(M ) (/ is compactly supported), and 

(1.5) for blow-up case, B~(a,/) : f(y) has a "minimal growth" as y — > oo, 

equations (II. 3p become two true nonlinear eigenvalue problems to study, which can be 
considered as a pair of mutually "adjoint" ones. All related aspects and notions used 
above, remaining still somehow unclear, will be properly discussed and specified. 
Our goal is to show by using any means that, for small n > 0, eigenvalue problems 



^1.6) (NEP)± admit countable sets of solutions $ ± (n) = {az, f~}\<y\>o, 



where 7 is a multiindex in M N to numerate the pairs. 
The last question to address is whether these sets 

(1.7) $ ± (n) of nonlinear eigenfunctions are evolutionary complete, 

i.e., describe all possible asymptotics as t — > +00 and t — > CT (on the corresponding 
compact subsets in the variable y in (II. 2p ) in the CP for the TFE-4 (11.11) with bounded 
compactly supported initial data. 
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Our main approach is the idea of a "homotopic deformation" of fll.l[) as n — > + and 
reducing it to the classic bi-harmonic equation 

(1.8) u t = -A 2 u in R N xR + , 

for which both problems (jl.6p and (11. 7p are solved positively by rather standard (but not 
self-adjoint) spectral theory of linear operators. 

1.2. Main TFE applications: nonnegative and oscillatory solutions. It has been 
well known since the 1980s, when higher-order parabolic models began to be studied 
more actively, that the TFEs-4 like (II. ip have many applications arisen particularly in 
modeling the spreading of a liquid film along a surface, where u stands for the height 
of the film in this context of thin film theory. Other physical related problems come 
from lubrication theory, nonlinear diffusion, flame and wave propagation (the Kuramoto- 
Sivashinsky equation and the extended Fisher-Kolmogorov equation), phase transition at 
critical Lifshitz points and bi-stable systems. We refer to a number of key survey and 
other papers on TFE theory such as [3j HJ [8j El [TOj [121 [30] ; see also Peletier-Troy [35] as 
a guide to higher-order ODEs and [HI [21] for most recent short surveys and long lists 
of references concerning physical derivations of various models, key mathematical results 
and applications of TFEs. Concerning mathematics of TFEs, one has to refer to the 
pioneering Bernis-Friedman paper [5] and [7J [23j H3] for the role of source-type similarity 
solutions of (II. ip . On modern existence-uniqueness theory for the ID TFE (for FBP 
setting), see [29], [201 § 6], and references therein. 

It should be pointed out that most of the results cited above are associated with non- 
negative solutions of a free-boundary problem (FBP) for the TFE-4 (II. ip . while currently 
this equation is written for solutions of changing sign. Moreover, as mentioned above, the 
development of general approaches to nonnegative solutions of the FBP began with the 
work of Bernis-Friedman [5] in 1990 with such solutions having a most relevant physical 
motivation and applications. 

The study of oscillatory solutions of changing sign for the TFE-4 is more recent; see 
[T3"l I2TJ1 [22] and references therein. It was shown in [IH] [22] (see also [T] as the most 
recent publication) that such solutions can be attributed to the Cauchy problem (CP) in 
R N x R + , rather than a FBP, posed in a bounded domain with moving free boundaries. 
The study of the Cauchy problem is interesting from both points of view in some biological 
applications as well as its clear mathematical interest in PDE theory. We refer to PQ, where 
more details on the CP setting are available. 

In this connection, another pioneering paper of Bernis-McLeod in 1991 [6] should be 
mentioned, where existence and uniqueness of first three oscillatory source-type solutions 
of the Cauchy problem for the fourth-order porous medium equation (PME-4) 

(1.9) u t = -{\u\ n u) xxxx in lxl + , 

are studied. Here, unlike (II. ip . equation (11. 9p contains a monotone operator in the metric 
of H~ 2 (R). By classic theory of monotone operators [33], the CP for (11.91) with compactly 
supported initial data uq admits a unique weak solution that is oscillatory close to the 
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interfaces for all n > and evidently for n — 0, where it becomes the bi-harmonic equation 
fll.Sp . with an oscillatory kernel of the fundamental solution; see below. 

For n > 0, such classes of the so-called "oscillatory solutions" of TFE-4 ( jl.ip is dif- 
ficult to tackle rigorously, and even their ODE representatives (in the radial geometry) 
exhibit several surprises in trying to describe sign-changing features close to interfaces, 
[19]. Indeed, the CP in ~R N x R + shows compactly supported blow-up patterns, which 
have infinitely many oscillations near the interfaces and exhibit maximal regularity there 
(consult [19j for further details). It turns out that, for a better understanding of such 
singularity oscillatory properties of solutions of (jl.ip , it is quite fruitful to consider the 
homotopic limit n — > + , thanks to the spectral theory developed for the pair {B, B*} in 
[T5] for rescaled operators where n = 0. Thus, here we perform a homotopic approach, 
more rigorous than before, in order to obtain such interplay between the CP for the TFE-4 
(jl.ip and the bi-harmonic equation (j 1.8ft . 

1.3. Our approach, problem setting, and layout of the paper. Before giving a 
description of our approaches, it is worth mentioning again that TFE theory for free 
boundary problems (FBPs) with nonnegative solutions is well understood nowadays (at 
least in ID). The FBP setting assumes posing three standard boundary conditions at 
the interface, and such a theory has been developed in many papers since 1990. The 
mathematical formalities and general setting of the CP is still not fully developed and 
a number of problems remain open. In fact, the concept of proper solutions of the CP 
is still partially obscure, and moreover it seems that any classic or standard notions of 
weak-mild-generalized-... solutions fail in the CP setting. 

Various ideas associated with extensions of smooth order-preserving semigroups are well 
known to be effective for second-order nonlinear parabolic PDEs, when such a construction 
is naturally supported by the maximum principle. The analysis of higher-order equations 
such as (jl.ip is much harder than the corresponding second-order equations or those 
in divergent form (cf. (11. 9ft ) because of the lack of the maximum principle, comparison, 
order-preserving, monotone, and potential properties of the quasilinear operators involved. 

It is clear that the CP for the bi-harmonic equation (jl.8p is well-posed and has a unique 
solution given by the convolution 

(1.10) u(x,t) = b(x --,t) * uo(-), 

where b(x,t) is the fundamental solution of the operator D t + A 2 . By the apparent 
connection between ( jl.ip and (I1.8P (when n — 0), intuitively at least, this analysis provides 
us with a way to understand the CP for the TFE-4 by using the fact that the proper 
solution of the CP for ( jl.ip . with the same initial data Uq, is that one which converges 
to the corresponding unique solution of the CP for ( 11 ,8p . as n — > 0. Thus, we shall use 
the patterns occurring for n = 0, as branching points of nonlinear eigenfunctions, so some 
extra detailed properties of this linear flow will be necessary. 

In Section [3j we, more carefully, introduce two classes of similarity solutions (the so- 
called nonlinear eigenfunctions), while Section H] is devoted to necessary properties of the 
spectral pair {B, B*} of linear differential operators that occur at n = 0. 
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Our further analysis is as follows: 



1.4. Local bifurcation-branching analysis for global solutions (Section |4]): first 
operator theory discussion. In the first part of this work, we perform a local bifurcation- 
branching analysis with respect to the continuation parameter n > 0, when that parameter 
is small enough. Thus, we obtain the bifurcation of solutions of the non-gradient equation 
(II. 3p , from the branch of the corresponding eigenfunctions of a rescaled linear operator. 
This yields some information and properties of the global in time similarity solutions 
dL2]) + of the TFE-4 (O]). 
The linear elliptic equation occurring at n = 0, 

(1.11) BF = -A*F + \yV y F + %F = in R N , J F(y) dy — 1, 

R N 

where F is the rescaled kernel of the fundamental solution b(x, t) in (ll.lOp . will be pivotal 
in the subsequent analysis. Indeed, the nonlinear operator in (jl.3p . 

(1.12) B+(a, /) := -V • (\f\ n VAf) + l=ft y .Vf + af, 
can be written in the following equivalent form: 

(1.13) B+(a, /) = -A 2 / + ±=p y ■ Vf + af + V • ((1 - |/| B )VA/). 

Then, the solutions of B+(oj, /) = are regarded as steady states of the nonlinear evolu- 
tion equation 

(1.14) / T = B+(a,/) in R N x R+. 

The bifurcation-branching point from our solutions (for n = 0) will be denoted by (n, /) = 
(O,"0fc), which is shown to occur for some values of the nonlinear eigenvalue a written as 
a = ak (k = |/3 1 is characterized by a multiindex /3 in R N ), and ipk representing the 
eigenfunctions of the operator B, whose expressions will be obtained in detail later on. 

Firstly, we shall prove that no bifurcation from the branch of trivial solutions (n, f) = 
(0, 0) occurs when the parameter n approximates 0. Secondly, an infinite number of 
branches of solutions is shown to emanate from the eigenfunctions of the rescaled linear 
operator B. Consequently, this analysis provides us with a countable family of solutions 
pairs (11.61) for the nonlinear equation (11. 3p , for small n > 0. 

According to classic bifurcation theory [151 EEl EH [38], we denote 

(1.15) B+(a, /) = F(n, f) := £(«, n)f + J\f(n, /), 

and assume that n is the main continuation parameter. Then, in order to have a branch 
of solutions emanating from the branch of trivial solutions (n, /) = (0, 0) at certain values 
of the parameter n (bifurcation points), the nonlinearity in (11.131) . denoted by 

JV(n,/):=V.((l-|/| B )VA/), 
must fulfill the following conditions: 

(1.16) (NL) : JV(n,0) = 0, D f M{n,0) = for all nGl + . 
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In other words, 

tffaf) = o(\\f\\) as f^O. 
Recall that, here, Af(n, f) serves as a perturbation of the operator B+ defined as in 
(11.131) . Thus, under the given assumptions (which are not that easy to pose in a suitable 
functional setting, to say nothing of the proof), the linear operator denoted by 

(1.17) £(a,n) := -A 2 + ^y V + al, 

defines an analytic semigroup in the space, where the solutions of (jl.3p , are defined. 

Note that, in any case, the necessary assumptions for the nonlinearity of (jl.3p , are far 
from clearly specified, when / is very close to zero, so something else must be imposed. Let 
us note that the condition (NL) in (I1.16p . roughly speaking, assumes that the functions 
f(y) are sufficiently smooth and have "transversal" zeros with a possible accumulating 
point at a finite interface only. Otherwise, if f(y) exhibits vanishing inside the support 
at a sufficiently "thick" nodal set, with many non-transversal zeros, this can undermine 
the validity of (11.161) . even in any weak sense. 

As customary in nonlinear operator theory, instead of the differential operators in (I1.15p . 
one has to deal with the equivalent integral equation 

(1.18) / = -(C(a, n) - aI)-\Af(n, f) + af), 

where a > is a parameter to be chosen so that the inverse operator (a resolvent value) 
is a compact one in a weighted space L 2 p {R N ); see Section [31 We will show therein that 
the spectrum of £ is always discrete and, actually, 

(1.19) a(C(a,n)) = {(1 - an) ( - f) + a, k = 0, 1, 2, ...}, 

so that any choice of a > such that a (jL cr(£) is suitable in fl!.18p . Therefore, in 
particular, the conditions (I1.16P are assumed to be valid in a weaker sense associated 
with the integral operator in ( 11. 18ft . 

Let us explain why a certain "transversality" of zeros of possible solutions f(y) is of 
key importance. As we see from (11.131) . we have to use the expansion for small n > 

(1.20) \f\ n - 1 = e n ln|/l - 1 = 1 + n In |/| + ... - 1 = n In |/| + ... , 

which is true pointwise on any set {|/| > Eq] for an arbitrarily small fixed constant Eq > 0. 
However, in a small neighbourhood of any zero of f(y), the expansion (11.201) is no longer 
true. Nevertheless, it remains true in a weak sense provided that this zero is sufficiently 
transversal in a natural sense, i.e., 

(1.21) J^i^ln|/| as n^0 + 

in L™ c , since then the singularity ln|/(?/)| is not more than logarithmic and, hence, is 
locally integrable in (II. 18ft . Equivalently we are dealing with the limit 

nln 2 |/|^0, as n 1 + , 

at least in a very weak sense, since by the expansion (11.201) we have that 

^i-ln|/| = lnln 2 |/| + .... 
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Note also that actually we deal, in (|1.18p . with an easier expansion 



(1.22) (|/r - 1)VA/ = (n In |/| + ...)VA/, 

so that even if f(y) does not vanish transversally at a zero surface, the extra multiplier 
VAf(y) in ( 11 .221) . which is supposed to vanish as well, helps to improve the corresponding 
weak convergence. Furthermore, it is seen from (I1.13P that, locally in space variables, the 
operator in f l 1 . 1 8 j) (with a = for simplicity) acts like a standard Hammerstein-Uryson 
compact integral operator with a sufficiently smooth kernel: 

(1.23) /-(VA^KI/r-^VA/]. 

Therefore, in order to justify our asymptotic branching analysis, one needs in fact to 
introduce such a functional setting and a class of solutions 

V = {f = f(;n): / 6 H 4 p (R N )}, 

for which: 

(1.24) V : (VA)- 1 (^iVA/) (VA^Qn |/|VA/) as n -» 0+ 

a.e. This is the precise statement on the regularity of possible solutions, which is necessary 
to perform our asymptotic branching analysis. In ID or in the radial geometry in ~R N , 
(ll.24p looks rather constructive. However, in general, for complicated solutions with 
unknown types of compact supports in WL N , functional settings that can guarantee fll.24p 
are not achievable still. We mention again that, in particular, our formal analysis aims 
to establish structures of difficult multiple zeros of the nonlinear eigenf unctions f>y(y), at 
which (ll.24p can be violated, but hopefully not in the a.e. sense. 

To study nonlinear integral operators it is necessary to construct a function space 
in which the integral operator possesses favorable properties (continuity, compactness). 
Indeed, one can apply the classical fixed point principles of Schauder's type to an operator 
acting between suitable Banach spaces. In this situation we can assert the existence of 
such fixed points establishing the continuity and boundedness of the integral operator. 
To do so, thanks to classical nonlinear integral operator theory we should impose the 
continuity of the kernel function involved in our integral operator (ll.23p . 

Within the previous context, let us observe that the integral equation with a Hammerstein- 
Uryson operator-type fl 1 . 2 3 [) is equivalent to the integral equation (I1.18p . for which we 
know that the inverse operator (C(a, n) — al)~ l is compact. Indeed, by the spectral theory 
described in Section [3j we are able to deduce that the operator C is defined between two 
exponential weighted spaces. Hence, it looks like to ascertain the existence of such fixed 
points for (11.181) and, equivalently for (11.231) . the suitable Banach spaces (that will pro- 
vide us with the existence of solutions of the original equation (I1.12p ) are precisely those 
exponential weighted spaces, together with the assumption of continuity of the kernels 
involved in the equivalent integral equations (I1.18P and (11.23)) . 

In addition, we would like to mention that for the study of elliptic problems of order 
2m by Schauder's inversion procedure the suitable Banach spaces could be the typical 
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pairs consisted of either Holder spaces 

(C 2 p m ' a {R N ),C° p a {R N )), with 0<a<l, 
or as remarked in the previous paragraph Sobolev spaces 

(Wp m ' p (R N ), L p p (R N )), with 1< p < oo. 

The particular weights assumed for those Banach spaces should be consistent with 
the exponential ones obtained in section [3j Thus, this enables us to obtain a priori 
estimates for the solutions of the original nonlinear equation ( I1.12p and provides us with 
the compactness of the integral operators involved in (11. IS)) and (ll.23p . 

1.5. Regularity convention. Overall, we observe that, unlike classic existence bifurcation- 
branching theory [161 EH EE], where sufficiently smooth expansions are used, the present 
singular one ( 11.20)) dictates a special functional setting in a subset V of functions (admis- 
sible solutions), for which (11.24)) should be valid a priori. In particular, such an analysis 
of the integral equation ( II . IS)) will always require some deep knowledge of admissible 
structure of proper solutions f(y) near zero (nodal) sets (also unknown), which we are 
still not aware of. Recall that, as our main goal, the present branching analysis is going 
to give us a first understanding of such delicate properties via the known eigenfunctions 
of the linear rescaled operators to appear at n = + . 

Thus, since these necessary nodal properties of possible solutions f(y) are unknown 
entirely rigorously, we perform our analysis under the following regularity convention: 

we assume that, regardless of the strong degeneracy of the nonlinear elliptic operator 
involved, the problems under consideration in both integral and differential forms 
admit sufficiently regular expansions of solutions in small n > 
in the functional class, for which (11.24)) holds. 

As usual in bifurcation theory, the hypothesis for this to be valid is formulated for 
the equivalent integral representation of the operators, though, for simplicity, we perform 
the n-expansion analysis in the simpler (but indeed equivalent) differential form. Overall, 
currently, we honestly do not think that our analysis can be justified more rigorously than 
that: technicalities to arise can be extreme and a full prove truly illusive. However, in 
Appendix A, we show that a suitable justification of the branching is indeed achievable 
provided that clear transversality of a. a. zeros of linear eigenfunctions of B is known. 
Nevertheless, the problem of the actual existence of nonlinear eigenfunctions for small 
n > remains open still. 

1.6. Further branching discussion. Now, once we have discussed the principal dif- 
ficulties, which have arisen, we carry out a local bifurcation analysis close to n = 0. 
Then, the change of stability from the branch of trivial solutions (n, /) = (n, 0) should 
be determined by the spectrum of the linearization C(a,n) and the assumptions of the 
nonlinearity. Therefore, a bifurcation would take place at some values of the parameter n 
if every neighbourhood of (n, /) = (n, 0) in R x C(Q) contains a nontrivial solution (n, f) 
of (1 1 . 3 p under the assumptions imposed for the linear and nonlinear part of the equation. 
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However, as will be proved later, such a bifurcation from the branch of trivial solutions 
(n, /) = (n, 0) never happens at the value of the parameter n — 0. Hence, we obtain a 
branching from the points (n, /) = (0, ipk) only; these arguments are fully consistent with 
more particular results obtained earlier in pp. 

Moreover, since the bifurcation from the branch of trivial solutions depends on the 
eigenvalues of the linear operator (11.171) . we believe that such a bifurcation does not exist 
at all (any proof is also very difficult). Indeed, after some rescaling of the type y i-> ay, 
with a = (1 — cm)~ 1//4 , it turns out that the spectrum of (11.171) is directly related to the 
spectrum of the linear operator B. This suggests that no bifurcation from the branch of 
trivial solutions ever happens. Therefore, as also discussed in p], we conjecture that if 
one wants to ascertain the global bifurcation analysis for these similarity TFEs-4, a new, 
different operator theory approach must be used. 

Next, let us obtain the values of the parameter a for which the bifurcation-branching 
phenomena occurs. The spectrum of the operator B in (II. lip , which appears for n = 
in (II. 3p , . is already well known and will be explained in detail in the next sections. This 
has the form 

a(B) = {A, =: -f , = 0,1,2,...}. 

Moreover, for k = 0, i.e., for the first eigenvalue-eigenfunction pair {a (n),f }, from 
the conservation of mass condition, denoting by M(t) the mass of the solutions of (II. ip . 
we have that (here Q is the rescaled support of f(y), however, for the CP, one can put 

n = R N ) 

M(t) := j u(x, t) dx = t~ a j /(f) dx = t~ a+pN j f(y) dy. 
n n a 

This yields the exact values 

(1.25) - a + (3N = =► a (n) = ^ and o (n) = 

However, the construction of the first eigenfunction fo{y) is not that straightforward 
even in ID; see [20, § 7], where its oscillatory properties cease to exist at a heteroclinic 
bifurcation calculated numerically as 

n h = 1.7587.... 

It is worth mentioning that, fortunately, for all n G (0, 1) (this interval is of particular 
interest in what follows), both the existence and the uniqueness of fo(y) follow from the 
results of [6], since, rather surprisingly, source-type similarity profiles for (II. ip and (II. 9p 
are reduced to each other with the parameter change n i— > ^-; see a precise statement in 
[201 Prop. 9.1]. 

Thus, it turns out that, when the parameter n approximates zero, we obtain according 

to (H25D 

«o(0) = f , 

so the solutions of (II. 3p , seem to approach the first eigenfunction ip associated with 
the first eigenvalue of the operator B, i.e., corresponding to Ao = 0. However, that 
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approximation for the solutions of (jl.3p , should also be extended to the eigenfunctions 
ipki for any k > 1, when the parameter a reaches the following values: 

(1.26) afc (0):=-A fc + f for any k — 1,2, ... , 

where \ k are the eigenvalues of the operator B, so that 

ao(O) = f , «i(0) = a 2 (0) = . . . , a k (0) = ^ . . . . 

Then, we introduce the next expression for the parameter a 

(1-27) a fc (n) := ^ - A fc . 

Hence, due to the necessary assumptions, the structure of the bifurcating-branching set 
emanating at (n, f) = (0,ipk) depends on the spectral theory for the operator (11. lip . 
For the first eigenvalue, since Ao = is simple, we can ascertain accurately the local 
bifurcation-branching. Then, at least for sufficiently small n's, the bifurcation-branching 
is locally a C 1 curve, which can be parameterized asm (n (s), f(s)) in R x C(Q) with 

Mo), /(o)) = (o, Vo), f(o) = $0, $o e n, 

where ' := 4-, emanating from the eigenfunction ip at the n = in the direction of the 
space Yq orthogonal to the eigenspace ker B, with x/jq — F (the rescaled fundamental 
kernel of b(x,t)) being the eigenfunction associated with the eigenvalue Ao = 0. 

However, in the case when the multiplicity of the eigenvalues A& with k > 1 is higher 
(bigger than 1), we obtain that the continua emanating at (n, f) = (0,tpk) are tangent 
to the manifolds Y k , orthogonal to ker (B + And, hence, we might have more 

than one direction of bifurcation-branching, depending on certain values related to the 
eigenfunctions which generate the eigenspace. This certainly agrees with the work of 
Rabinowitz [36], in which for potential operators and bifurcation from the branch of trivial 
solutions, one of the next alternatives for the bifurcation structure must be obtained: 

(i) for the value of the parameter where the bifurcation takes place, the trivial solution 
is not isolated; or 

(ii) for any other value of the parameter in one-sided neighbourhood of the bifurcation 
point, there are at least two nontrivial solutions; or 

(iii) for any other value of the parameter in a neighbourhood of the bifurcation point, 
at least one nontrivial solution exists. 

In general, the question about how many precise branches bifurcates for any k > 1 
remains an open problem, though we think it is very related to the dimensions of the 
eigenspaces. As far as we know, only partial and very specific results have been obtained 
for non-variational problems with higher multiplicities. 

1.7. Blow-up patterns via branching theory (Section H]). This is a natural coun- 
terpart of the global similarity analysis of PDEs in the limits as t — > oo. We next consider 
blow-up limits as t — > T~ < oo, or t — > 0~ as in (11. 2p . where T = 0. We thus perform 
a detailed and systematic analysis of the blow-up similarity solutions. This is done again 
by using the homotopic approach as n — > + via branching theory, but this time based 
on the Lyapunov-Schmidt methods in order to obtain relevant results and properties for 

10 



the solutions of the self-similar equation (jl.3p . This homotopic-like approach is based 
upon the spectral properties of the adjoint (to the B above) operator 

(1.28) B* = -A 2 -±yV, with a(B*) = {A^ = -f, |/3| = 0,1,2,...}, 

which occurs after blow-up scaling of the linear counterpart (jl.8p of the TFE-4 (11.11) 
for n — 0. Note that (11.281) admits a complete and closed set of eigenfunctions being 
generalized Hermite polynomials, which exhibit finite oscillatory properties. 

It is curious that, in [19], blow-up similarity analysis of the related unstable TFE-4 did 
not detect any stable oscillatory behaviour of solutions near the interfaces of the radially 
symmetric associated equation. All the blow-up patterns turned out to be nonnegative, 
which is a specific feature of the PDE under consideration therein. This does not mean 
that blow-up similarity solutions of the CP do not change sign near the interfaces or 
inside the support. Actually, it was pointed out that local sign-preserving property could 
be attributed only to the blow-up ODE and not to the whole PDE (11.11) . Hence, the 
possibility of having oscillatory solutions cannot be ruled out for every case. Indeed, 
thanks to the polynomial expressions of the eigenfunctions for the operator B*, we have, 
in particular, that the first eigenfunction ip*(y) = 1 is not oscillatory, but for some other 
eigenfunctions we shall expect sign-changing behaviour. 

Then, this homotopy study exhibits a typical difficulty concerning the desired structure 
of the transversal zeros of solutions, at least for small n > 0. Proving such a transversality 
zero property is still a difficult open problem, though qualitatively, this was rather well 
understood in ID and radial geometry, [T9] . 

1.8. TFE: FBP and CP problem settings. We recall that, for both the FBP and the 
CP of (II. ip . the solutions are assumed to satisfy standard free-boundary conditions: 

{u = 0, zero-height, 
Vw = 0, zero contact angle, 

— n ■ V(|ii| n Ait) = 0, conservation of mass (zero-flux) 

at the singularity surface (interface) T [u} = dQ, which is the lateral boundary of 

supp u C R N x R + , N>1, 

where n stands for the unit outward normal to r [u], which is assumed to be sufficiently 
smooth (the treatment of such hypotheses is not any goal of this paper). For smooth 
interfaces, the condition on the flux can be read as 

lim n- V(\u\ n Au) = 0. 
dist(x,r M)4-0 

For the FBP, dealing with nonnegative solutions, this setting is assumed to define a 
unique solution. However, this uniqueness result is known in ID only; see [29] , where the 
interface equation was included into the problem setting. We also refer to [201 § 6.2], where 
a "local" uniqueness is explained via von Mises transformation, which fixes the interface 
point. For more difficult, non-radial geometries in M , there is no hope of getting any 
uniqueness for the FBP, in view of possible very complicated shapes of supports leading to 

n 



various "self-focusing" singularities of interfaces at some points, which can dramatically 
change the required regularity of solutions. 

For the CP, the assumption on nonnegativity is got rid of, and solutions become oscil- 
latory close to interfaces. It is then key that the solutions are expected to be "smoother" 
at the interface than those for the FBP, i.e., f l 1 . 2 9 [) are not sufficient to define their reg- 
ularity. These maximal regularity issues for the CP, leading to oscillatory solutions, are 
under scrutiny in [20]; see also [Tj, as the most recent source of such a study. 

Next, denote by 

M(t) := / u(x,t) dx 

the mass of the solution, where integration is performed over smooth support (1R is al- 
lowed for the CP only). Then, differentiating M(t) with respect to t and applying the 
divergence theorem (under natural regularity assumptions on solutions and free bound- 
ary), we have that 

r n{t} 

The mass is conserved if J(t) = 0, which is assured by the flux condition in (ll.29p . 

The problem is completed with bounded, smooth, integrable, compactly supported 
initial data 

(1.30) u{x, 0) = u (x) in F [u] n {t = 0}. 

In the CP for ( 11. lft in M. N x R + , one needs to pose bounded compactly supported initial 
data f l 1 . 3 f) prescribed in WL N . Then, under the same zero flux condition at finite inter- 
faces (to be established separately), the mass is preserved; however smoother regularity 
properties of solutions require a separate study/understanding; see [20] for some results. 

2. Self-similar solutions: two nonlinear eigenvalue problems 

2.1. Global similarity solutions. We now more carefully derive the problem for global 
self-similar solutions of (II. ip . which occur due to its natural scaling-invariant nature. 
Namely, using the following scaling in (11. ip : 

x := fix, t := At, u := uu, with 

(2.1) 

du v du du v du d u v d u 

dt A dt ' dxi ii dSi ' dx'f /i 2 dx'f ' 

and substituting those expressions in (II. ip yields 

if = -^v-(i«rvA^). 

To keep this equation invariant, the following equalities must be fulfilled: 

_,n + l _ x R , 4/3-1 

(2.2) 



1 = ^ /i := ^> v := An, so that 

u(x, t) := X~^~ u(x, t) = X~^~ «(-), where t — A. 



Consequently, we have to rescale in the following way: 

X 



(2.3) u + (x,t) = t~ a v(y,r), y := f , r = In t : E + — )■ R, where f3 = l=f±, 
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such that f(y, r) = u(3r, r), we obtain, after substituting (I2.3P into (ll.ip and rearranging 
terms, that / solves a quasilinear evolution equation given by 

(2.4) v T = B+(a,v) = -V -(\v\ n VAv) + ±=fty-Vv + av in R N x R + . 

Consider the steady-states of the parabolic equation (12.41) . Thus, we analyze the local 
bifurcation-branching behaviour of the nonlinear eigenvalue problem: 



(2.5) B+(a,/) = -V.(|/rVA/) + ^y.V/ + a/ = 0, f e C (R N ) 



Here, the "boundary conditions at infinity" stated as / £ Cq(R. n ) are naturally associated 
with the known properties of finite propagation for TFEs, which have been mathematically 
justified about two decades ago at least; see a survey on energy methods in PDE theory 
in [28]. Then, any assumption stating that f(y) is "sufficiently small" at infinity, e.g., 

(2.6) feH A (R N ) or H*(R N ) 

(the last space is a domain of the linear operator B in ( 11.11)) ; see the next section), 
would lead to compactly supported solutions. In fact, such a conclusion entirely depends 
on asymptotic (i.e., local, not any global) properties of the nonlinear elliptic operators 
involved, so will not be a main concern in our study. 

2.2. Blow-up similarity solutions. The blow-up similarity solutions of the TFE-4 ( II. II) 
correspond to completely different limits and then describe a "micro-scale" structure of 
its solutions at any given point. For convenience, we reduce the blow-up time to T = 0. 
Then, similar to the global solutions, replacing 1 1-> (—t), we obtain the patterns 

(2.7) U-(x,t) := (—t)~ a f(y), y = jzkp, with the same parameter j3 = 1 ~ 4 n " . 

Hence, substituting that expression into ( 11. ip and rearranging terms, we arrive at the 
following quasilinear elliptic equation: 



(2.8) B-(a,/) = -V.(|/|*VA/)-A/.V/-c*/ = in 



In order to get the corresponding second "adjoint" nonlinear eigenvalue problem, one 
needs to specify a "minimal growth" of admissible nonlinear eigenfunctions f(y) as y — > 
oo. This will be done in Section |5] Note that, by obvious and straightforward reasons, 
( I2.8p does not admit compactly supported solutions. Indeed, the nature of blow-up scaling 
( 12 .7p would then mean disappearance of a such a solution in finite time, contradicting 
uniqueness and other easy asymptotic issues for the TFE-4. 

In general, for solutions with finite blow-up time T £ R, the full self-similar scaling 

(2.9) u(x, t) = (T- t)- a w(y, r), y:= r = - ln(T - t) : (-oo, T) R, 
yields the parabolic equation 

(2.10) w T = B-(a,w) = -V ■ (M n VAw) - ^f±y ■ Vw - aw in l N xl + . 
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3. Spectral properties of the linear operator B 



In this section, we describe the spectrum cr(B) of the linear operator B obtained from the 
rescaling of the bi-harmonic equation (II. 8p . This spectral theory will be essentially used 
in ascertaining the direction of the branches bifurcating from the trivial (actually nonex- 
istence) and other eigenfunctions and the number of branches for the blow-up solutions. 

3.1. Relation to a linear eigenvalue problem. Let u(x,t) be the unique solution of 
the CP for the linear parabolic bi-harmonic equation (11. 8p with the initial data 

(3.1) u eL 2 p (R N ), where p(y) = e a ^ 4/3 , a > small, 
given by the convolution Poisson-type integral 

(3.2) u(x,t) = b(x,t) *u = t^T J F((x — z)t~*)u (z) dz. 

Here, by scaling invariance of the problem, the unique fundamental solution of the oper- 
ator + A 2 has the self-similar structure 

at 



(3.3) b(x,t)=t--F(y), y:=-^ (x E 



Substituting b(x,t) into (11.81) . we obtain that the rescaled fundamental kernel F in 
(13. 3p solves the linear elliptic problem (II. lip . B is a non-symmetric linear operator, 
which is bounded from H p (W N ) to L 2 p (W N ) with the exponential weight given in (13. ip . 
Here, more precisely, a G (0, 2d) is any positive constant, depending on the parameter 
d > characterizing the exponential decay of the rescaled kernel: 

(3.4) \F(y)\ < De~ d M 4/3 in R N (D > 0), 

Later on, by F we denote the oscillatory rescaled kernel as the only solution of (II. lip , 
which has exponential decay, oscillates as \y\ — > oo, and satisfies the standard pointwise 
estimate (13 .4p . 

Thus, we need to solve the corresponding linear eigenvalue problem: 



(3.5) B^ = Xip in R N , ijj e L 2 (R N ). 



It seems clear that the nonlinear problem (II .ip formally reduces to (13. 5ft at n = with 
the following shifting of the corresponding eigenvalues: 



(3.6) A = -a(0) + 



N 
4 ' 



It is another reason to call (12. 5 j) a nonlinear eigenvalue problem, since for n = it reduces 
to the classic eigenvalue one for a linear differential operator. Moreover, crucially, the 
discreteness of the real spectrum of the linear problem (13. 5p can be apparently inherited 
by the nonlinear one, but a complete justification of this issue is far from being clear. 
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3.2. Functional setting and semigroup expansion. Thus, we solve (I3.5P and calcu- 
late the spectrum of cr(B) in the weighted space L 2 (M N ). We then need the following 
Hilbert space: 

H*(R N ) c L 2 p (R N ) c L 2 {M N ). 
The Hilbert space Hp(M ) has the following inner product: 

4 

(v,w) := / p{y) D k v(y)D k w(y)dy, 
where D k v stands for the vector {D@v , \(3\ = k}, and the norm 



\H 2 P := J p(y) E \D k v(y)\ 2 dy. 

Next, introducing the rescaled variables 

(3.7) u(x,t)=t-T W (y,T), y:=&, r = kit : R + R, 
we find that the rescaled solution w satisfies the evolution equation 

(3.8) w T = Bw, 

since, substituting the representation of u(x, t) (13. 7p into (jl.8p yields 

-A> + \ y ■ V y w + f w (y, r) = iff. 

Thus, to keep this invariant it must be satisfied that = 1 ==>- r = Int. Hence, w(y, r) 
is the solution of the Cauchy problem for the equation (13.81) and with the following initial 
condition at r = 0, i.e., at t — 1: 

(3.9) w (y) = u(y, 1) = 6(1) * u = F * u . 

Thus, the linear operator J^— B is a rescaled version of the standard parabolic one J^ + A 2 . 
Therefore, the corresponding semigroup e Bl ~ admits an explicit integral representation. 
This helps to establish some properties of the operator B and describes other evolution 
features of the linear flow. From (13.21) . we find the following explicit representation of the 
semigroup: 

(3.10) w(y,r) = f F(y — ze~5) u (z) dz = e Br w , where x = t*y, r = lnt. 

Subsequently, consider Taylor's power series of the analytic kernel 

(3.H) F(y-ze-i) = E W e-^ { -^D^F(y)z^ = E W e^^(y)^ 

for any y G Mr, where 

_/3 y Pi ~Pn 

and ipp are the normalized eigenfunctions of the operator B. The series in (13.1 ip converges 
uniformly on compact subsets in z G M N . Indeed, estimating coefficients for \/3\ — I, 



T^^DPFiy)^...^ 
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< b,\z\ 



by Stirling's formula we have that, for I ^> 1, 

(3.12) b t = ^sup yeRNm=l \D^F(y)\ « fl~ l / 4 e 1 / 4 « /"^V = e - /ln3/ / 4 + nnc . 

Note that the series J^&iM ^ as its radius of convergence R = oo. 
Thus, we obtain the following representation of the solution: 

(3.13) w{y, r) = £ (/3) e^M^Mv)^ wh ere A^ =: -f 

and {ipp} are the eigenvalues and eigenfunctions of the operator B, respectively, and 

R N 

are the corresponding moments of the initial datum wq defined by (13. 9p . 

3.3. Main spectral properties of the pair {B, B*}. Thus, the next results hold [T8] : 

Theorem 3.1. (i) The spectrum o/B comprises real eigenvalues only with the form, 

(3.14) *(B) := {A, =: -Ifl , \fi\ = 0, 1, 2, ...}. 

Eigenvalues \p have finite multiplicity with eigenfunctions, 

(3.15) := t^^Ffe) s (&)* ' ' • m- 



(ii) T/ie subset of eigenfunctions $ = {^3} is complete in L 2 (WL N ) and in L 2r ' ,: " A ■ 

(iii) For any A ^ f(B), the resolvent (B — A/) -1 is a compact operator in L 2 (R N ). 

Then, the adjoint operator B* of B (in the dual metric of L 2 (M. N ) takes the form HI . 28[) 
and is defined in the weighted space L 2 *(R N ), with the domain H%(R N ), where the (dual) 
weight function is exponentially decaying: 

p*^ = m = e ~ alyr > °- 

It is a bounded linear operator |18j . 

B* : Hp,(R N ) ->• L 2 , (R N ), so (Bv,w) = (v,B*w) , v G H}(R N ), w G H}.(R N ). 

Moreover, the following theorem establishes the spectral properties of the adjoint operator 
which will be very similar to those shown in Theorem l3.1l for the operator B. 

Theorem 3.2. (i) The spectrum o/B* consists of eigenvalues of finite multiplicity, 

(3.16) o(B*) = o(B) := {A,, =: -f , = 0, 1, 2, ...}, 

and the eigenfunctions ip%(y) are polynomials of order \(3\. 

(ii) The subset of eigenfunctions $* = is complete and closed L 2 ,(R N ). 

(iii) For any A ^ cr(B*) the resolvent (B* — A/) -1 is a compact operator in L 2 *(1Sl N ). 
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It should be pointed out that, since ijj = F, 

f RN ipo dy = f RN F(y) dy = 1. 

However, thanks to (13 . 1 5[) we have that 

J r n^/3 = for any ^ 0. 

This expresses the orthogonality property to the adjoint eigenfunctions in terms of the 
dual inner product. Due to Theorem l3.2l the adjoint eigenfunctions are polynomials which 
form a complete subset in L?*(M. N ) with exponential decaying weight p*(y) = e _a ' y ' 4/3 . 

Note that [IS], for the eigenfunctions {^V?} of B denoted by (13.1 5p . the corresponding 
adjoint eigenfunctions are generalized Hermite polynomials of the form 

[0/4] 

(3-17) ^(V):=^[V /J + £ J! A V]- 

Hence, the orthogonality condition holds: 

(3.18) (i/jp,^) = 5p„ for any /3, 7, 

where (•, •) is the duality product in L 2 {R N ) and <5^ i7 is the Kronecker's delta. Operators 
B and B* have zero Morse index (no eigenvalues with positive real parts are available). 
The main spectral results are extended [18] to 2mth-order linear poly-harmonic flows 

(3.19) u t = —(—A) m u in R N x R + , 

where the elliptic equation for the rescaled kernel F(y) takes the form 

(3.20) BF = -(-A y ) m F+^yV y F+£ l F = in R N , f F(y) dy — 1. 

In particular, if m = 1 and N = 1, we find the classic second-order Hermite operator B 
(see [TT] for further information) 

BF = F" + I F'y + i F = 0, 

whose name is associated with the work of Charles Hermite of 1870, although such equa- 
tions and polynomial eigenfunctions of the adjoint operator B* = Dy — | yD y were ob- 
tained earlier by Jacques C.F. Sturm in 1836, [57] : see [2U Ch. 1] for history and references. 

4. (NEP) + : LOCAL BIFURCATION-BRANCHING ANALYSIS VIA A FORMAL APPROACH 

In this section, we show the nonexistence of local bifurcations points from the trivial 
solution (n, f) = (n,0) for the nonlinear operator (j!.13p . and, hence, the existence of 
branching from the eigenfunctions of the linear operator B when n is sufficiently close to 
zero. This analysis allows us to show locally the existence of non-zero solutions of (|1.3p | 
in the proximity of n — 0. 

Throughout this section, we write the operator (I1.13P in the form Fin, f), denoted by 
(I1.15p . where C(a, n) and J\f{n, f) are the corresponding linear and nonlinear parts of the 
operator (I1.13P under the abstract framework already explained above (first section). This 
operator is of class C r , with r sufficiently big to make all the subsequent derivatives exist. 
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Moreover, when n = we have the operator B defined by (11.111) . for which we showed in 
the previous section its complete spectral theory. It is apparent that its eigenvalues 
and eigenf unctions ipk, with k > 0, will determine the precise number of branches from 
(n,/) = (0,V*). 

Let us note that the nonlinearity condition assumes that the functions / are sufficiently 
smooth and have "transversal" zeros with a possible accumulating point at a finite inter- 
face only. Also, for each k > 0, denoting uq± = 0, we find that 

(4.1) ker[£(n 0i A;)] = ker (B + f /) = span {ifrp, \(3\ = k} for any k = 0, 1, 2, 3, • • • , 
where £(n ,fc) := £(atk{0), no,*)- Then, due to Fredholm's alternative (see e.g., [16]). 

R[C(n 0j k)} = {u E C(Cl) : f n uipk = 0}, such that 
ker (B + f/) ®R[C(n 0tk )] = L 2 p {R N ) for any k = 0, 1, 2, • • • . 

Thus, ipk R[C(n 0ik )] for any k > 0. It is clear that the operator £(n 0jfc ) = B + 1 1 is 
Fredholm, i.e., i?[£(a,n)] is a closed subspace of L^R^) and 

dimker(£(a, n)) < oo, codim_R[£(a, n)] < oo, 

at least for each n m + . Then, the operators £(rio,fc) are Fredholm of index zero. Indeed, 
we already know that the first eigenvalue Ao = is a simple one of the operator £(n ,o) = 
B, so its algebraic multiplicity is 1. Hence, we will apply the classical results of Crandall- 
Rabinowitz [15] about bifurcation for simple eigenvalues in order to prove the nonexistence 
of bifurcation points from the branch of trivial solutions at the value Uq = 0. 

Note that, due to Theorems 13.11 and I3.2[ for any k > 1, the algebraic multiplicity is 
equal to the geometric ones, so we are not dealing with the problem of introducing the 
generalized eigenfunctions (no Jordan blocks are necessary for restrictions to eigenspaces) . 

On the other hand, when dealing with essentially non-analytic functions of n, at n = 
as in (I1.13p . we cannot use standard apparatus of bifurcation-branching theory even in 
the case of finite regularity; cf. [TU [311 [38] . This reflects the main partially technical but 
often principal difficulties of such a branching study. 

Once the assumptions are established, we introduce the following concept, which will 
play a role in the forthcoming analysis. In dynamical system theory, such concepts are 
typical for characterizing various types of bifurcations. 

Definition 4.1. (n ^, 0), with n 0j fc = for any k = 0, 1, 2 • • • , is a bifurcation point for 
equation (12. 5p from the curve of trivial solutions (n, 0), if there exists a sequence 

(n km J m )eRx(H*(R N )\{0}), 

where m > 1, such that 

lim m ^.oo(n fem ,/ m ) = (0,0) and J r (n fcm ,/ m ) = for each m > 1. 
Since B + 1 1 is Fredholm, and 

(4.2) D f T{n Q)ki Q)f=(B + il)f 
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for any k, on the whole, as was discussed in [31], it is clear that the condition 

M fc =ker (B + fl) > 1, 

is necessary for any (n, f) = (^o,fc,0) to be a bifurcation point of (12. 5 j) from (n, 0), the 
trivial solution. However, the bifurcation does not occur depending only on the linear 
part. Therefore, as was shown in [31] through simple algebraic examples, the nature of 
the nonlinearity determines the sufficiency condition in order to have such a bifurcation. 
To be more precise, the nonlinear part of the operator (I1.13P must fulfill the assumptions 
established in the first section of this paper. 

Therefore, according to bifurcation theory, the points, where a bifurcation from the 
branch of trivial solutions occurs, must depend on the spectrum associated with the linear 
part and assuming some conditions on the nonlinear part. Here, we prove the nonexistence 
of bifurcation from the branch of trivial solutions at the value of the parameter n = 0. 
Moreover, due to the rescaled relation between the linear operators B (11.1 ip and C(a,n) 
defined by (I1.17p . it turns out that there is no bifurcation point from the branch of trivial 
solutions at any n > 0. 

4.1. Bifurcation-branching for simple eigenvalues. Firstly, we present a result that 
provides us with the nonexistence of the branch emanating from the trivial solutions at 
the point {n, f) = (n o ,fc,0), with n ok = 0, when X k is a simple eigenvalue. Secondly, 
consistent with some recent findings [TJ, we also show that there exists a branching from 
the eigenfunction ip at the value of the parameter no,o = 0. We actually know that Ao = 
is a simple eigenvalue of B in a general setting in Mr, but we cannot assure that it is 
the only such one in other geometries. For instance, in ID and in the radial setting, all 
eigenvalues are simple, so we can apply this simplified analysis. 

Thus, the calculus below will be valid for any k such that A& is simple (under suitable 
restrictions). However, to avoid excessive notation, we make all the computations for the 
case k = 0, which is always special and simpler. 



Lemma 4.1. Under the regularity convention and assumptions in Section [Lj 



(i) (n, /) = (rio,o, 0) ; with no,o = ; is not a bifurcation point for the stationary equation 

F{n, f) = (/ e C (R N ) or H*(R N ); cf (Q). 

(ii) (n, f) = (no,Oj V'o) ^ s a branching point for the stationary solutions of the functional 
J-{n, f). Furthermore, letY be a subspace of Hp(R N ) , 

Y := {u E C(R N ) : J RN uip = 0}, such that ker(£ ,o) © *o = H*{R N ). 

Then, there exists e > and two maps of the class C r ~ x , 

n Q : (-£,£) ->• R, $ : {-e,e) -»• Y , n (0) = 0, $ (0) = 0, 

such that, for n Q = and for each s G (— e, e), 



(4.3) 



F{no{s), f (s)) = 0, f (s) := + s$ (s), 
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where ipo is the eigenf unction associated with the simple eigenvalue Ao of B, so that 
dimker[£(no)] = 1. Moreover, there exists p > such that ifj r (n,f) = and (n, /) G 
-B r (0, ?/>o), i/ien either f = ip Q , or (n, f) = (n (s), fo(s)) for some s G (— s,e), where 
B r (0,i/)o) is the a ball of the radius r centered as (0,ipo) in IR 2 x L 2 (M. N ). Furthermore, if 
T is analytic, so aren (s), ao(s), and fo(s) nearO. 

Proof (i) To show that (n, /) = (n , 0), with n ,o = 0, is not a bifurcation point, we check 
that the transversality condition of Crandall-Rabinowitz [15] is not satisfied. Under the 
regularity assumptions and convention of Section [L4l imposed on the linear and nonlinear 
parts, we obtain that C{n) := DfFin, 0), where standard calculations of the derivative 
are allowed: 

D f F{n, f)g:= lim^ HnJ +h g)~HnJ) 

_ i- -V-(\f+h9\ n VA(f+hg))+i^y-V(f+hg)+ a (f+hg)-T(nJ) 

— nm h^0 ^ 

_ ,. T(nJ)+h[-Vin\f\ n ' 1 gVAf)~Vi\f\ n VAg)+^y-Vg+ag]+o(h)-T(n,f) 

— Uni/i^o ji 

= -V ■ (n|/| n -^VA/) - V ■ (|/ | n VA<?) + l=f±y.Vg + ag, 

for any g G Cq(M. n ), or Hp(M N ); cf. (12.61) . Moreover, owing to the spectral theory shown 
in Section El we find that there exists a singular value Ao, the eigenvalue of the operator 
£( n o,o)> with n fi = 0, associated with the eigenfunction ip . Set £ ,o := £(^o,o) an d 
£i, : = £^( n o,o) where 

d_r(n\ — ( - N(i+Nn)-N 2 n „ _ N 2 j 

dn^VV ' — I 4(4+Vn) 2 h 4 V (4+7Vn) 2 J 1 

Then, ker(£ ,o) = spanj-^o} and the following transversality condition does not hold: 
(4.4) A,c*Mfl[£o,o]. 



Indeed, suppose 



£isfll>o = -Jr V^ - C G i?[£ ,o] , so 



(4.5) _g I/ . V ^ ) -g^ = -A 2 t; + Jy.Vt; + f «, 

for any v G i?[£ ,o]- Hence, multiplying (14. 5 p by the adjoint eigenfunction i[)q and inte- 
grating by parts yields 



16 J ruru 1 16 J ^uw - ro 1 16 



which implies the nonexistence of bifurcation at no,o = from the trivial solution. 
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(4.6) G(s,n ,$( 



(ii) We now prove the final statement of Lemma 14.11 Then, under the same necessary 
regularity assumptions and the convention in Section IL~4l we define the auxiliary operator 

.r(no,i/>o+s$o) ^ / q 

s ' ' ' 

D f F(n , ^ )$o, if s = 0, 

for s G R, s ~ 0, no e R, and $o G ^o- Since J 7 is C r in both variables, Q is C r_1 in all 
its arguments. The number r is sufficiently large ensuring that the derivatives employed 
in the sequel exist. Moreover, by the definition, we have that 

(4.7) £(0,0,0) = 0, 

since D/J r (no, 0)$o — by construction. Then, if the zeros of the eigenfunction ipo are 
transversal a.eQ, we find that 

1-M,P W «0 

in the weak sense (or even "a.e.") for a sufficiently small s. Hence, 

D M) g(o, o, o)K $o) = hm^o 8(w y (w) 

(4.8) = lim- - Rl^( n oM(h^ ) 



h->a jt 

lim^ (£ ,o + /m £i,o)$o + o(h) = £ ,o$o- 



Thus, owing to 
the operator 



y ©ker(£ 0i o)=# p 4 (R 



D Mo) g(Q,0,0) : HxY^Uj^ 



is an isomorphism. Then applying the implicit function theorem, we deduce the existence 
and uniqueness of two C r_1 functions 

n : (—£,£)—»• R, $ : (—£,£)—>■ Y , such that 
7i (0) = 0, $ (0) = 0, Q{s, n (s), $„(*)). 

□ 

Throughout the rest of this section, we calculate the possible types of local bifurcations. 
According to our formal analysis above, it follows that (12. 5 j) has a local curve of solutions 

(n (s), /o(s)), /o(s) : = V>o + s$ (-s), 

emanating from (n, /) = (n, t/>q) at n = 0. This shows in some natural sense that the 
component emanating at (n, /) = (0, ip ) consists of two subcontinua, £q and <£q in the 
direction of $ an d — $0; respectively, where $ belongs to the orthogonal space Y to the 



1 For ipaiy) which is a radial function, this is very probable; however we do not have a fully convincing 
proof. 
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eigenspace spanned by the eigenfunction associated with the simple eigenvalue A = 0. 
Moreover, those functions (no(s), fo( s )) admit the next expansions of the form: as s — > 0, 

(A9) ^o(s) :=S7o,i + s 2 7o, 2 + o(s 2 ), 

/o(s) := + s$ ,i + s 2 %, 2 + o(s 2 ), 

for certain real numbers 70,/ and some functions $o,« £ ^O; with I = 1,2. In addition, since 
we are assuming that a is dependent on n, this eventually yields 

"o( s ) = f + sr/0,1 + s 2 ?7o,2 + o(s 2 ), 

where 770,2 6 ^ f° r ari Y i = 1,2, ■ ■ ■ . 

Now, substituting the expansions defined by (14. 9 p and the corresponding expansion for 
the parameter ao(n) into the nonlinear elliptic equation (12. 5p yields 

-V- (|Vo + s$o,i| (s7o ' 1+o(s)) VA^o) -sV- (|Vo + ^o,il (s7o ' 1+ ° (s)) VA$ ,i) 
_ S 2 V . () ^ o + s $ 0>1 |(*>m-k>(')) vA$ 0i2 ) + \y ■ V(^ + s$ 0) i + 

-s— 4 — y ■ V(^o + s$ ,i + o{s)) 

+ (f + sr/0,1 + s V2 + o(s 2 )) (V»o + s$ ,i + o(s)) = . 

Next, passing to the limit as s — >• 0, by the regularity convention (in particular, this 
assumes the transversality condition for the zeros of the eigenfunction ^ ), we have that 
|^ | s to,i =i_|- S j Q1 l n |-0 O | _|_ o(s), and according to the spectral theory shown in Section[3] 
we find that 

B^j = (- A 2 + ±y V + f /)^o = 0. 
Hence, dividing the rest of the terms by s and passing to the limit as s — > gives 

(4.10) B$ ,i = 7o,i [g y V^o + V • (ln |^ | VA^o)] + r/0,1^0 • 

Now, applying Fredholm's theory [16J to ( I4.10p yields that there exists a function <&o,i> 
which solves (14.101) if and only if the right hand side is orthogonal to kerB, i.e., to the 
eigenfunction = 1 of the adjoint operator B*. This uniquely defines the coefficient: 

(AH) 'V ■= VQ,i(n^o) = rjo^ 

l ' ' m ' 1 (V *,S/-V^o) + (^,V-(ln|^ |VA^o)) ^(l, J /-VVo>+(l,V-(ln|^o|VAVo)> ' 

provided that the denominator does not vanish (notably, a difficult property to prove or 
even to verify numerically). Therefore, any different branching-type in the vicinity of 
n = from the first eigenfunction ip Q of the operator B will depend on the values of the 
coefficients 7 0j i and 770,1 related by (14.111) . 

4.2. Bifurcation-branching for semisimple eigenvalues. Hereafter, in this section, 
we focus on the case when the kernel is multidimensional. Note that, the question of 
local bifurcation at the value of the parameter, for which the corresponding eigenvalue is 
simple has been extensively studied in literature. In particular, as was shown above, for 
(12. 5p . under some assumptions over the nonlinearity, we proved that no bifurcation takes 
place at n = from the branch of the trivial solution, and when the eigenvalue is simple. 
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However, for any neighbourhood around n = 0, a branch of solutions emanates from the 
associated eigenfunction in the direction of the orthogonal subspace Y . 

On the other hand, for eigenvalues with higher multiplicity, we prove that the nontrivial 
solutions emanating from the eigenf unctions ip k at the value of the parameter no,fc = 0, 
for any k > 1, are tangent to a manifold Y k . 

Also, in general, it is not completely understood how many branches emanate from the 
trivial solution, which remains an open problem and it can only be obtained for some 
specific examples. 

For the case of bifurcation from the branch of trivial solutions (n, f) = (n, 0), there 
exist some results supposing that the operators are potential (see [T7] and [3S]) and very 
few for non-gradient, non-self-adjoint operators, [52] . 

Here, we provide a number of possible branches of bifurcation-branching when n is close 
to + from the eigenf unctions ip k under some conditions imposed over certain values. 

Similarly to the case of simple eigenvalues, we already know that 

(4.12) M k = dimker (B + \ I) > 1, 

for any k > 1, and the inequality is strict in dimensions N > 1. Note that A^ = — | is 
not a simple eigenvalue and (14. ip is fulfilled by ip k as the eigenfunctions associated with 
those semisimple eigenvalues of B, A/., such that ip k '■= J2\p\=k cp^P-, f° r every k > 1 and 
under the natural "normalizing" constraint 

(4-13) £1,1=^ = 1. 

Subsequently, under the circumstances imposed for the nonlinearity, it is apparent that 

dimker (B + \l) = codimi?(B + 1 1) = M k , 

for any k > 1. Thus, as was discussed above, the operator (14. 2 p is Fredholm of index zero 
since it is a compact perturbation of the identity. 

It should be pointed out that the odd crossing number condition might fail, so we are 
not distinguishing between odd or even multiplicities (see the works Ambrosetti [2] for 
gradient operators and by Kromer-Healey-Kielhofer [32J for more general operators). It 
is classically known that, when the multiplicity is odd, there is always a bifurcation. How- 
ever, when the multiplicity is even, the bifurcation depends strongly on the nonlinearity. 

The next theorem is one of the main results of this paper. 

Theorem 4.1. Let the assumptions for the linear and nonlinear part of the functional J 7 
be satisfied, together with the regularity convention in Section fTm and (I4.12p hold. Then: 



(i) ( n o,k = 0, 0), for any k > 0, is not a bifurcation point, and; 

(ii) if ' 

ker (B + *l)®Y k = H}(R N ), 
where the subspace Y k is defined by 

Y k := {ueC(R N ) : J RN u^J k = 0}, 
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there exists e > and two maps of class C r 1 , 

n k : (-e,e) ->■ R, : (-e,e) V*, 

swc/i £/ia£, for n = and /or eaca s G (— e, e), 
(4.14) ^Ms),/*(s)) = 0, :=V* + s$a(s). 

Furthermore, if J-{n, f) is analytic in a neighbourhood of(0,ipk)> k — 1,2, so are n k (s), 
a k (s), and fk(s) near s = 0, a countable number of branches emanate from ip k forn « + . 

Note that our earlier result for simple eigenvalues, A; = 0, is included here, with similar 
conclusions. 

Proof. Firstly, we consider the following auxiliary operator: 



(4-15) g(s,n k ^ k ):-- , . 

for s G R and close to zero, G R and G Y^. Since J 7 is C r in both variables, Q is 
C r_1 in all its arguments. As customary, we impose some regularity conditions making 
sure that all the derivatives in the sequel exist. Moreover, by definition we have that 

(4.16) £(0,0,0) = 0, 

for every k > 1. Thus, similarly as done for the case of simple eigenvalues (I4.8p . 

D (nM g(0, 0, 0)K, $ fe ) = hm^ g(o.n fc ^)-g(o,o,o) 
(417) = hm^ 2lZ^t^ 

= \im h ^ (C 0tk + hn k Ci yk )<&k + o{h) 

= £>o,k®k- 

Hence, if the following condition 



n©ker(£o, fc )=^ p 4 ( TOA, 



holds, then the operator D( nfcj $ fc )C/(0, 0, 0) : R x Y k — >■ L^(R ) is an isomorphism. Con- 
sequently, we can apply the implicit function theorem. Therefore, the existence and 
uniqueness of the following two C r ~ l functions are guaranteed: 

n k : (— e, e) — > R, $fc : (s,s) — > Y k , such that 

Now, in order to conclude the proof, we must show that if (n, f) = {n o k ,0), with 
n o,k = 0, is not a bifurcation point, for any k > 0, then the following condition, providing 
us with the bifurcation from the branch of trivial solutions at the value of the parameter 
no k = 0, must not be satisfied 



span ,£i,fc$Aij k }©-R(B + \l) =L 
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where Ci^ := 4->C(no,&), with no,k = 0, and {^i, • • • , ipM k } a basis of the subspace ker (B+ 
\ /) such that 

= Ci^i H h C Mk 1pM k , 

with the "normalizing" constraint (14. 13ft . Thus, we suppose that 

A,*^* = cjCxjdj = (-f + ^f)y ■ Eft cjVipj - § Eft Cj^j e #[£o,fc], so 
f4 . Eft q,A A(= -S + * )v • Eft ^ - § Eft c,i- 

for any t> 6 i?[£o,fc]- Now, we restrict ourselves to the case when k — 1 and = 2 
(i.e., iV = 2) to avoid excessive calculations. Hence, multiplying (I4.18j) by the associated 
adjoint eigenfunctions "0* and ip 2 and integrating by parts, we obtain the following system: 

^±1 ffcy ■ ( Cl V^i + c 2 V</> 2 ) + § / ^( Cl ^ + c 2 ^ 2 ) = 0, 

= 0, 



/ 4*y • ( Cl V& + c 2 V^ 2 ) + § / ^(c^x + c 2 ^ 2 ) 



and, hence, 



(4.19) 



jV+l 
16 

jV+l 
16 



+ c 2 
+ c 2 



Af+l 
16 

Af+l 
16 



= 0, 
0. 



Consequently, if the determinant of the system (14.191) for the unknowns C\ and c 2 is 
different from zero, 



^0, 



(N+l)j ■ V^i + N 2 J fcfa (N + 1)J r x y ■ V^ 2 ) + N 2 j foj> 2 
(N + l)J i,*y • + iV 2 / ^ (JV + 1) JT ^ ■ V^ 2 + iV 2 / 4^ 2 

we arrive at the desired result with the normalizing constraint (I4.13p . Similar computa- 
tions can be done for any finite k and M^; see below. □ 

Furthermore, as was performed for the case with simple eigenvalues, we ascertain the 
conditions that provide us with how the branching from the eigenfunctions at n = is 
and how many branches we actually have. Hence, by Theorem 14. 11 for any k > 1, the 
local curve of solutions 

(njfe(s), fk(s)), fk(s) ■= i>k + s$fc(s), 

emanates from the branch of solutions (n, f) = (0,^), for any k > 1. That curve of 
solutions is defined locally by two maps of class C r_1 , 

n k : (-£,£) ->• E, : (-£, e) -> Y" fc , 

such that, for s = 0, 

n fe (0)=0, $ fc (0) = 0, 

and the eigenfunction of the subspace ker (B + | /) , as well as the expansion of the 
parameter a(n) depending on n. 
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Moreover, since the eigenvalues associated with the eigenfunctions ^ are semisimple, 
the dimension of the kernel M k will be greater than 1. Thus, the component emanating 
at n = will do it in the direction of the orthogonal manifold to the one generated by 
• • ■ ,ipM k }, in such a way that ip k = c\ip\ + • ■ ■ + CM h ^M k - In other words, depending 
on the coefficients ci, • • • ,CM k , we shall obtain different directions of the bifurcation- 
branching. Then, those functions (rifc(s), fk(s)) admit the following expansions of the 
form: as s — > 0, 

n k (s) := sj kjl + s 2 7 fci2 + o(s 2 ), 



fk{s) := Ipk + S$ fc ,i + S 2 $ fej2 + 0(S 



2\ 

) 



for certain real numbers 7^ and some functions ^ k ,i £ Yk, with / = 1,2, and k > 1. 
Furthermore, we set ctfc(s) = + s^7a,-,i + s 2 ?7fc,2 + °( s2 ) ? where 77^ G 1R for any z = 1, 2, • • • 
and any k > 0. 

Hence, substituting those expansions into the equation (12. 5ft and dividing by s gives 

-V • (|^ + s$ M | (s7fc ' 1+o(s)) VAV> fc ) - sV ■ (|^ + s$ M | (s7fc ' 1+o(s)) VAd> M ) 
_ S 2 V . ()s ^ + s 2 $fcji |( S 7M+°W) V A$ fc , 2 ) + | r V(^ + + o(s)) 

-s A 2 '- — y ■ v(ipk + s$ k ,i + o(s)) 

+ (™ + S7? M + s V 2 + o(s 2 )) {i) k + s$ M + o(s)) = . 

Then, passing to the limit as s — > in a similar way as was done for the case of simple 
eigenvalues (assuming a "sufficient transversality" of a. a. zeros of the eigenfunctions ipk 
and, hence, the expansion \ip k \ slk - 1 = 1 + £7^1 In \ipk\ + we have that 

(B + |/)^= (- A 2 + ±^ + ^7)^ = 0, 

which is true by Section [3j Dividing the rest by s and letting s — > give 

(4.20) (B + I J)$ M = 7M y • V^ fc + V • (In \^ k \ VA^)] + , 

where ^ = Ci^i + • • • + CM k ^pM k - Hence, multiplying by the associated adjoint eigen- 
functions l^*, • • • , V'm. }, integrating over WL N and applying the Fredholm alternative [16], 
we obtain an algebraic system with the coefficients {cj, j = 1, • • • , M k }, 7^,1, and rj^i as 
the unknowns. Again, to avoid excessive calculations, we will restrict ourselves to the 
simplest case in which k = 1 and Mi = 2 (N = 2). Then we arrive at the following 
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algebraic system: 

ci P 1 ' 1 ^ /^t 1/ ■ Wi + 71,1 / ^* V • (In \cA + c 2 i> 2 \ VA^O 

+7i,i M* V • (In |ci^i + c 2 4| VA^ 2 ) + t?!,! / ^* 4] = 0, 
( , .21 ) <; Cl [7M(g±fc) j ^ y . v ^ + 7i)i f ^ V • (ln | Cl ^ + C2 ^ 2 | VA^x) 

+7/1,1 / fcti] +02 / & y ■ W 2 

+7i,i Ji>2 V • (In |ci^i + c 2 4| VA^ 2 ) + t?!,! / ^\ $ 2 ] = 0, 
Cl + c 2 = 1. 

We will achieve the existence of solutions due to standard fixed point theory arguments 
(see [1] for further details). Then, in order to ascertain how many possible solutions we 
might have, one can fix the value of 77^,1 and solve the nonlinear algebraic system (14.21 j) for 
the remaining unknowns. Thus, from the third equation of (I4.2ip we find that C\ = 1 — c 2 . 
Then, setting ci^i + c 2 -?/> 2 = ipi + c 2 (t/>2 — ^1), substituting it in the other two equations, 
and integrating by parts in some of the terms (the ones with the logarithm), we obtain 

7i,i [t ffcvWi+J Vfa h] + 771,1 / ^ 

+c 2 [2iig±fi ffcy. ( V^ 2 - V$i) + 71,1 / (h 2 -h 1 ) + 771,1 / r x - fa)] = 0, 
7i,i [t Ikv + / W 2 * /ij + 771,1 / 4* & 

+c 2 [^ /^y(V^-V^i)+ 7 i,i JV^I^-^+r?!,! M*(^-^i)] =0, 

where /i x := In \fa + c 2 ('02 - wl VA^i , h 2 := In \tp 1 + c 2 (^ 2 - t/>i)| VA-^ 2 . 
Hence, we next solve the system without considering the extra perturbation terms, which 
h\ and h 2 are involved in, i.e., 

w<(ti,i,C2) := 71,1 / V^J" /ii + Ca7i,i J V^* (h 2 - h) , with z = l,2. 

Hence, we need to solve the system, 

C2 [^^1 J fay. ( - Wi) + 771,1 / & (4 - &)] 

(4 22) +7i,i^ M v ■ v& + VIA I k k = 0, 

C2 [IMg±a /&y.(V&-Wi)+ifc,i 

+71,1^ / 4* y • v^i + 771,1 / 4* fa = 0. 

At this point it is quite easy to prove that, for example, after substituting the expression 
for c 2 , obtained from the first equation, into the second equation, we arrive at a quadratic 
form which depends on the unknown 71,1, 

5(7i,i) := ^i7i 2 ,i + 5i 7 i,i + Ci = 0, 

with at most two possible solutions. Note we are assuming that at least one of the 
following conditions is fulfilled: 

(4.23) IhimV ffty. (yfa - V^i) + r/i,! / fa { (fa -fa)^0 with i = 1, 2. 
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Those solutions will correspond to two possible values of c 2 which are the roots of the 
following quadratic form: 



0(ca) := A 2 c\ + B 2 c 2 + C 2 = 0. 

Moreover, owing to the "normalizing" constraint ( I4.13p . we have that c 2 £ [0, 1]. Hence, 
for that quadratic form the following is ascertained: 

(i) C2 = 05(0) = C 2 ; 

(ii) c 2 = 1 (5(1) = A 2 + B 2 + C 2 ; 

(iii) Differentiating (5 with respect to c 2 , we obtain that &'(c 2 ) = 2c 2 A 2 + B 2 . Then, 
the critical point of the function <S is c\ = — ^ and ©(c^) = — ^ + C 2 . 

Once the solutions for c 2 are established, that we know they are between and 1 according 
to the "normalizing" constraint, we are able to ascertain the solutions for j 2 ,i- Although, 
they can reach any value in the real line, and not only values between and 1, they 
must fulfill the quadratic form $ as well, so that, we will have at most two solutions 
corresponding to this unknown 72,1. Therefore, going back to ( 14 . 2 3 [) and supposing it is 
true, we shall obtain two solutions after imposing the following conditions: 

(a) C(A 2 + B 2 + C 2 ) > 0; 

( b ) C(--fc + C 2 ) <0; and 

(c) < -it 2 < 1- 

On the other hand, if — j^- + C 2 = then we have just one solution of the quadratic form. 

However, in the case when condition (I4.23P is not fulfilled, we obtain a single unique 
solution for j 2 ,i, which satisfies the following equality: 

_ via Jj|ji _ 771,1 J i>i h. 

Unfortunately, in this case nothing can be said about the number of solutions for c 2 and, 
hence, for the other unknowns appearing in the system (14.211) . unless just one of them is 
satisfied, in which case we ascertain one unique solution for all the unknowns. Observe 
that, for any solution pair (71,1, c 2 ), there corresponds a value of 771,1, that we fixed above. 

Therefore, we will obtain at most two solutions of the system (I4.22p . and, eventually, 
imposing some conditions on the extra nonlinear terms £^(7^1, c 2 ) such as 

11^(71,1,02)11^ < min^^), 6(4)}, for any i = l,2, 

where, Ti 1 > c 2 are the values where the quadratic forms catch the critical points, we finally 
obtain at most two solutions for the original nonlinear algebraic system (I4.2ip . 

For the sake of completion, we extend these results to the case in which k = 2 and the 
dimension of the kernel is M 2 = 3 (again, N = 2). In other words, the kernel will be 
generated by {ipl, ip 2l V^} such that ip 2 = c\ipl + c 2 ip 2 + 03^3. Hence, for this particular 
case, multiplying again (I4.20p by the adjoint eigenf unctions ?/>*, with i = 1,2,3, and 
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imposing the "normalizing" constraint (j4.13|) . the following algebraic system is obtained: 
ci [ 72 ' l( X 6 +fc) / 4 y ■ V4 + 72,1 ffo V • (In IdV-! + c 2 4 + c 3 4|VA4) 

+%i / 4 4] +c 2 p*g±a / 4 y • V4 

+72,1 / 4 V ■ (In |ci4 + c 2 4 + c 3 4|VA4) + J/2,1 / 4 4] 

+72,1 / 4 V ■ (In | Cl 4 + c 2 4 + c 3 4| VA^) + %,i / 4 4] = 0, 

ci [^±^ ffcy V4 + 72,1 Ik V • (In |c x 4 + c 2 4 + c 3 4| VA^) 

+r/ 2ll /4*4]+c 2 p^ ikwk 

u „ ,x +72,1 / 4 V " (In |d^i + C2-02 + C 3 ^ 3 | VA^ 2 ) + J/2,1 / 4 4] 

1 j +c3p^M*</-v4 

+72,i / 4* V • (In | Cl 4 + c 2 4 + c 3 4| VA4) + / 4* 4] = 0, 

ci [ 72 ' l( X 6 +fc) M* 2/ • + 72,1 / 4 V • (In |c x 4 + c 2 4 + c 3 4| VA^) 

+w/4*4]+c 2 [^ ± ^i'4y-v4 

+72,1 / 4 V " ( ln l c i4 + c 2 4 + c 3 4|VA^ 2 ) +772,1 /44] 

+72,1 / 4 V ■ (In |ci$i + c 2 4 + c 3 4|VA4) + / 4 4] = 0, 
Ci + c 2 + c 3 = 1. 

As mentioned above, for the case fc = 1, the existence of non-degenerate solutions is 
guaranteed by standard fixed point theory. Moreover, since by the third equation C\ = 
1 — c 2 — c 3 , substituting it into the other three equations of the system f)4.24p yields 

72,1 Ikv vk + / v4 hi] + r)2,i I k k 

+C2 [IM^lS J^* y . (V4 - V&) + 72,1 / Vk (to - h) 

I k (k - k)] +c 3 p+mv Ikv (vk - v4) 
+72,1 / vk (h - h) + n2,i I k (k - k)] = 0, 

72,1 / 4* 2/ • V4 + / V4 ^l] + 7/2,1 / 4 4 

+C2 [z^g±a j 4* y . (v^ - vk) + 72,1 / vk (tkt - h) 
+r/ 2 ,i / 4 (4 - k)] +c 3 [^7^ / k y ■ (v4 - v4) 
+72,1 / v4 (to - to) + 772,! / k (k - k)] = 0, 

72,1 / 4* 2/ • V4 + / Vk hi] + 772,1 / 4* k 

+C2 phig±*i jfcy. ( v4 - v4) + 72,1 / vk (h 2 - to) 
+V2,i I k (4 - 4)] +c 3 Iky (vk - vk) 

+72,1 / vk (to - h) + m,! I k (4 - 4)] = 0, 
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(4.25) 



where, 



hi := In | -01 + c 2 (?/>2 - 0i) + c 3 (?/>3 - 0i)|VA-0i , 
/i 2 := In |^i + c 2 (^ 2 - i'l) + c 3 (^ 3 - ^i)|VA^ 2 , 
h 3 := In | -01 + c 2 (^ 2 - 0i) + c 3 (^3 - 0i)|VA-0 3 . 

Subsequently, as previously done for the particular case k — 1, we solve the nonlinear 
algebraic system ( I4.25P without including complicated nonlinear perturbations, which the 
terms hi, h 2 , and h 3 are involved in: 

w f (72,i s c 2 , Cs) := 72,1 / Vt/>* /ii + c 2 7 2 ,i / W* (h 2 - h) + c 3 7 2jl / V^* (/i 3 - fci) , 
with i = 1,2. Thus, ascertaining the number of possible solutions for the system 

72,1 MyV^+772,! ffcfa 

+c 2 [i^mhi jfcy. ( V ^ 2 - V0O + r/2,1 / 01 (02 - 0l)] 

+c 3 ffoy. (v4 - V&) + r/ 2il / # (4 - &)] = o, 

+c 2 [^lg±*i / 4* y • (V& - V&) + r7 2il / & - ^l)] 
+c 3 [^if±^ J 0* y • (V0 3 - V0O + 77 2>1 / 4* (03 - X )] = 0, 

72,1^ Ji>* 3 yVi>i+V2,i J Mi 

+C 2 F&!g±*l / 03* y • ( V^ 2 - V0!) + %1 / & (0 2 - &)] 
+C 3 / 4* y ■ (V^ - V&) + 77 2jl / r 3 - &)] = 0, 

and controlling the oscillations of the extra nonlinear perturbations cjj(7 2> ic 2 , C3), for any 
i = 1, 2, 3, as above for = 1, we achieve the desired results imposing the conditions 

/0* y .(V0 j -V0 1 )+77 ljl /^(^--fe^O, with 7 = 1,2,3 and j = 2,3. 

Our system can be reduced to the study of two perturbed quadratic forms. Therefore, 
we arrive at the problem of studying the number of intersections of two conic surfaces, 
which provides us with the number of solutions between zero and four. We postpone 
explaining how this approach works until Section 15.51 where it is applied to the blow-up 
nonlinear eigenvalue problem (jl.3p . This approach is quite similar for both the cases. 

As a preliminary but a key conclusion, it is worth mentioning now that we believe that, 
since we are dealing with a kernel of the dimension 3, in this case, we have four solutions. 
It seems then that two of them should coincide. This is a principal issue, owing to the fact 
that, somehow, the number of solutions depends on the coefficients we have for the system 
and, at the same time, on the eigenfunctions that generate the subspace ker (B + | i) , 
so its dimension. However, a full justification is not proved here and, due to the difficult 
nature of the problem, perhaps it will never be possible to justify it completely. 
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4.3. A short discussion on global behaviour of n-branches. After performing a 
very precise analysis about the bifurcation-branching analysis in the proximity of n = 0, 
we intend to explain how the global behaviour of the branch of solutions that emanates 
from the eigenfunction ipo at the value of the parameter n = can be. It is well known 
that the existence of solutions is not guaranteed for any n. Indeed, it was discussed in 
[20] that the existence of oscillatory solutions ends when n = 1.7587... . The value where 
the existence of oscillatory solutions ceases obviously depends on the type of the thin 
film equation we are dealing with (the pure TFE, with extra absorption terms, stable, 
unstable, etc). 

Due to the analysis performed in this section, we already know that there is no bifur- 
cation from the branch of trivial solutions at n = for the TFE (II. 3p , . 

B+(a,/) = -V-(|/| n VA/)+/32/.V/ + a/ = in R N . 

Moreover, after some rescaling the spectrum of the linear counterpart of (11.31) , is directly 
related with the spectrum of the operator B denoted by (II. lip 

BF = -A 2 y F + \y- V^F + f F = in R N , J F(y) dy — 1. 

R N 

Therefore, according to this, we believe that there is no bifurcation from the branch of 




> n 



Figure 1 . A formal global branching n-diagram. 

trivial solutions at any value of the parameter n > 0. However, a rigorous proof of this 
fact is difficult, since our nonlinear operators are not monotone and, hence, the main 
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techniques usually used in the analysis of second-order operators are not applicable here. 
However, if that was the situation, we could obtain the existence of certain turning points, 
or even a connection with another branch among the ones emanating from some other 
eigenfunction ipk, with k > 1; cf. Figure HJ We plan to carry out this work in a subsequent 
paper. 

5. Blow-up similarity profiles for the Cauchy problem via tvbranching 

5.1. Preliminaries: homotopy and nodal sets. In this section, we describe the be- 
haviour of the blow-up similarity solutions ( II. 2p of the TFE-4 (II. ip through the same 
homotopic approach by setting n \. in (jl.3p and, hence, arriving at the linear adjoint op- 
erator ( ll.28p . Then, we shall use the eigenfunction patterns occurring for n = (those are 
generalized Hermite polynomials ( I3.17P ) as branching points for nonlinear eigenfunctions 
providing us with a straightforward and practical rz-continuity approach to the self-similar 
equation ( 12. 8 p associated with the TFE-4 ( II. ip from the equation ( II. lip associated with 
the bi-harmonic equation (II. 8p . 

It is worth recalling now that homotopic approaches are well-known in the theory 
of vector fields and nonlinear operator theory (see [161 EI] f° r details). In our case, a 
"homotopic path" just declares the existence of a continuous connection (a curve) of 
some nonlinear eigenfunctions / = f k {y) satisfying ( 12 .8p that ends up at n = + at the 
linear adjoint polynomials ipl(y) given in (I3.17p . Due to Section [31 we already know that 
those profiles correspond to generalized Hermite polynomials given by ( I3.17p . which have 
finite oscillatory properties. For instance, for any even \/3\, the polynomials ( I3.17P do 
not have any zero nodal surface at all. However, for k > 1, linear combinations of such 
eigenfunctions do have nodal sets of known and relatively simple structure. 

For odd \/3\ (or on mult i- dimensional eigenspaces) , arbitrary linear combinations of 
Hermite polynomials for such a fixed k — \/3\ > 1 explain all possible structures of nodal 
sets and (see ES for a full formulation) 



(5.1) multiple zeros of solutions of the bi-harmonic equation (11. 8p . 



Furthermore, it turns out that, using classical branching theory, "nonlinear eigenfunc- 
tions" fk{y) of changing sign, which satisfies the nonlinear eigenvalue problem ( 12. 8p (with 
an extra "radiation-minimal-like" condition at infinity to be specified shortly), at least, for 
sufficiently small n > 0, can be connected with the adjoint polynomials ipl(y) in (I3.17p . or 
their linear combinations from the eigenspace. We are capable of justify this through the 
corresponding Lyapunov-Schmidt branching equation, trying to be as rigorous as possible 
in supporting and deriving the critical nonlinear eigenvalues «fc(n). 

5.2. Towards "minimal growth at infinity". This is about the "minimal" (a "radiation- 
like") condition at infinity, announced in (11.51) . which makes the equation (12. 8p to be a 
nonlinear eigenvalue problem. We recall that, for n = 0, the (NEP)_ in (11.31) reduces to 
the linear one for the operator B* in (ll.28p . with the straightforward correspondence 

(5.2) a k (0) = X k = -|, k — 0, 1,2, ... . 
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Equation (I2.8P admits two kinds of asymptotics at infinity. The first one is nonlinear 
and is given by the first two operators: assuming simple radial behaviour (/ = y 1 ) yields 



(5.3) -V.(|/rVA/)-i^j/.V/ + ... = f(y)~\yfi as y -> oo. 

Note that, as n — > + , precisely this behaviour leads to an exponentially growing bundle, 
which is prohibited in Theorem 13.21 by specifying the proper weighted space L^(M. N ) and 
eventually leading to the polynomial eigenfunctions (13.17j) . 
The second asymptotics is linear: as y — > oo, 

(5.4) -kfHy.Vf- a f + ... = f(y)~\y\\ where 7 = > 0, 

since by (15. 2p we have to assume that a^n) < (the first eigenvalue ao(n) = is not of 
particular interest; see below) and always dk{n) < -. Note that then 

(5-5) 7 - < i 

so that the linear behaviour ( 15. 4p is the actual minimal one in comparison with ( 15. 3ft . 

Overall, this allows to formulate such a "radiation-like" condition at infinity, which now 
takes a clear "minimal nature" : 



(5.6) find solutions f(y) of ( 12. 8ft bounded at infinity by functions as in ( 15. 4ft . 



In self-similar approaches and ODE theory, such conditions are known to define simi- 
larity solutions of the second kind, a term, which was introduced by Ya.B. Zel'dovich in 
1956 [39J. Many of such ODE problems (but indeed, easier) have been rigorously solved 
since then. For quasilinear elliptic equations such as ( 12. 8ft . the condition ( 15. 6 p is more 
subtle and delicate indeed. We cannot somehow rigorously justify that the problem (12. 8p . 
( 15. 6p is well posed and admits a countable family of solutions and nonlinear eigenvalues 
{ap(ri)}. We recall that using the homotopy deformation as n — > + was our original 
intention in order to avoid such a difficult "direct" mathematical attack of this nonlinear 
blow-up eigenvalue problem. 

We begin our actual study by noting that the first nonlinear pair for ( 12. 8p . (15. 6p is 
trivial: for any n > 0, 

(5.7) a (n) = and f (y) = 1, 

so that this well corresponds to the first Hermite polynomial from (13.171) with \/3\ = 0, 
where ipoiv) = 1- However, similarity solutions ( 12. 7p with the first eigenfunction in (15.71) 
are trivial and do not change sign, so, to understand formation of nonlinear "multiple 
zeros", we will study branching of eigenfunctions f^(y) for k > 1. 

5.3. Technical bifurcation calculus. Thus, the critical values ak{n) are obtained for 
small n > according to spectral theory established in Section [3J As was noticed, the 
explicit expression for the eigenvalues and eigenfunctions for the operator B* in (I1.28P 
are known; see Theorem 13.21 Moreover, supposing the corresponding linear counterpart 
from (13.171) with n = 0, we find, at least formally, that 

(5.8) C(a)f := -A 2 / - \ y ■ V/ - af = 0. 
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This equation can be considered as a linear perturbation in terms of the parameter a 
of that for the adjoint operator B* in (11.281) . From that equation combined with the 
eigenvalue expressions obtained for the operator B*, we derive the critical values for the 
parameter a given in (15. 2ft . where A& are the eigenvalues defined in Theorem l3.2[ Note 
that those eigenvalues coincide with the eigenvalues of the operator B. In particular, 
when k = 0, we have that «o = = Ao and the eigenfunction is i/jq = 1, satisfying 

B*^q = 0, so that ker£(a ) — spanj-^Q = 1}. 

Hence, Ao = is a simple eigenvalue for the operator £(azo) = B* and its algebraic 
multiplicity is 1. In general, we find that (note that k = is trivial) 

(5.9) ker (B* + f /) = span {ip^, \/3\ = k}, for any k = 0, 1, 2, 3, . . . , 

where the operator B* + | / is Fredholm of index zero. In other words, R[C(ak)} is a 
closed subspace of L 2 p (M. N ) and 

dimker(£(a fe )) < oo, codimi?[£(a fc )] < oo 

for each a k - Moreover, dimker (B + | J) = Mjt > 1, for any k — 0, 1, 2, 3, 

Then, once the relation between (15.81) and the linear operator B* has been established, 
for which we know its spectral theory, by regularity issues in Section 11.41 we can assume 
for small n > in (I3.17P the following expansions: 

(5.10) a k (n) := a k + ^ k n + o(n), \f\ n = \f\ n = e nln ^ := 1 + nln |/| + o(n), 

where the last one is assumed to be understood in a weak sense. Again, it is convenient to 
discuss further the last one. Indeed, the second expansion cannot be interpreted pointwise 
for oscillatory changing sign solutions f{y), though now these functions are assumed to 
have finite number of zero surfaces (as the generalized Hermite polynomials for n = 
do). However, as usual, this, of course, imposes some restrictions on the possible zeros of 
the eigenfunctions ip%(y)- According to the spectral theory in Section [3] we already know 
that those eigenfunctions and their linear combinations for the adjoint operator B* are 
generalized Hermite polynomials given by (I3.17p . Hence, they are analytic functions with 
isolated zeros. 

Since the possible zeros are isolated, they can be localized in arbitrarily small neigh- 
bourhoods. Indeed, it is clear that when |/| > 5 > for any 5 > 0, there is no problem 
in approximating of |/| n as in (15.101) . i.e., \f\ n = 0(n) as n — > + . However, when |/| < 5 
for 5 > sufficiently small, the proof of such an approximation is far from clear unless 
the zeros of the /'s are all transversal in a natural sense. In view of the expected finite 
oscillatory nature of solutions f~(y), this should allow one to obtain a weak convergence 
as in (11.211) to be used in the integral equation similar to (jl.ISp (with B replaced by B*). 

However, let us stress again that, in the present "blow-up" case, we do not need such 
subtle oscillatory properties of solutions close to interfaces, which are not known in com- 
plicated geometries. The point is that, due to the condition (15. 6p . we are looking for 
solutions f(y) exhibiting finite oscillatory and sign changing properties, which are similar 
to those for linear combinations of Hermite polynomials (13.171) . Hence, we can suppose 
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that their zeros (zero surfaces) are transversal a.e., so we find that, for n > and any 
S — 5(n) > sufficiently small, 

n|ln|/||»l, if \f\<S(n), 

and, hence, on such subsets, f(y) must be exponentially small: 

|ln|/||»I =► In |/| « -I => |/|«e-i 

Thus, we can control the singular coefficients in (I5.10p . and, in particular, see that 

(5.11) In |/| 6 LUR N ). 

Recall that this happens also in exponentially small neighbourhoods of the transversal 
zeros. 

It is worth recalling again that our computations below are to be understood as those 
dealing with the equivalent integral equation similar to (11.181) and operators, so, in par- 
ticular, we can use the powerful facts on compactness of the resolvent (B — A/) -1 and the 
adjoint one (B* — A/) -1 in the corresponding weighted L 2 -spaces. 

Note that, in such an equivalent integral representation, the singular term in (15.101) 
satisfying (15. lip makes no principal difficulty, so the last expansion in (I5.10P makes rather 
usual sense for applying standard nonlinear operator theory. Overall, the above analysis 
somehow justifies further branching study. We must admit that this is not a rigorous one, 
but is indeed sufficient for our formal expansions as n — > + . 

Thus, substituting (I5.10p into the nonlinear eigenvalue problem (12. 8p and omitting o(n) 
terms when necessary, we obtain the following expression: 

-V ■ [(1 + n In |/|) VA/] - ^V 1 ^ V ■ W - («* + A*i,*n)/ = , 
for any k — 0, 1, 2, 3, . . . . Hence, rearranging terms yields 

-A 2 / - nV ■ (In |/| VA/) - \y ■ V/ + ^±f^ y . V/ - a k f - ^ k nf = . 
In addition, using the expression of the operator B* yields 

(5.12) (B* + | /)/ + nN k {f) + o(n) = , 
with the operator 

(5.13) Af k (f) := -V • (In |/| VA/) + f r V/- » 1>k f . 

Subsequently, we shall compute the coefficients involved in the expansions ( 15. 10[) applying 
the classical Lyapunov-Schmidt method to (I5.12p (branching approach when n 10), and, 
hence, describing the behaviour of the blow-up solutions for at least small values of the 
parameter n > 0. Two cases are distinguished. The first one in which the eigenvalue 
is simple and the second for which the eigenvalues are semisimple. Note that due to 
Theorems 13.11 and 13.21 for any k > 0, the algebraic multiplicities are equal to the geometric 
ones, so we do not deal with the problem of introducing the generalized eigenfunctions 
(no Jordan blocks are necessary for restrictions to eigenspaces) . 
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5.4. Simple eigenvalue. Recall that this always happens for k = (not interesting) 
and also in ID and radial geometry, when all the eigenvalues of such ordinary differential 
operators are simple. 

As a typical example, we perform the analysis as for k = 0, bearing in mind the above 
other more interesting applications. 

Thus, since the first eigenvalue Ao = of B* is simple, the dimension of the eigenspace is 
Mq = 1, the analysis of this particular case presents less difficulties than the corresponding 
ones for any other k > 1. Hence, denoting ker B* = spanj-^Q = 1} and by Y Q * the 
complementary invariant subspace, orthogonal to if) , we set 

(5.14) f = r + v *, 

where V * G Y Q *. We define P * and P* such that P * + P* = I, to be the projections onto 
ker B* and Yq respectively. We next set 

(5.15) V* :=n${ + o(n). 

Then, after substituting (I5.14p into (I5.12p and passing to the limit as n — » + , we arrive 
at a linear inhomogeneous equation for Q\ 

(5.16) B*<^o = -MoM, 

since B*^q = 0. By Fredholm's theory [16] (spectral theory of the pair {B,B*} from 
Section [3] does also matter) , a unique solution V Q * G Y * of (15.16P exists if and only if the 
right-hand side is orthogonal to the one dimensional kernel of the adjoint operator, in this 
case B. In other words, in the topology of the dual space L 2 , the following holds: 

(5.17) (Af (r o )^o)=0. 

Therefore, ( I5.16P has a unique solution G Yq determining by (I5.15P a bifurcation 
branch for small n > 0. In addition, we obtain the following explicit expression for the 
coefficient /ii of the corresponding nonlinear eigenvalue ao(n) denoted by (15.101) : 

:= { ^ lRl tZ7 M0) = < V ■ (ln|^o|VA^),^ ). 

5.5. Multiple eigenvalues for k > 1. For any k > 1, we know that 

dim ker (B* + f l) = M* k > 1 (actually, M* k = M k ). 
Hence, we have to take the representation 
(5-18) f = E m=k c^ + V k \ 

for every k > 1. Currently, for convenience, we denote {^%}\p\=k — {^ii ■■•■> ^*M k }i ^ ne 
natural basis of the M k - dimensional eigenspace ker (B* + | J) and set ipl = ^1^1=* c /3^- 
Moreover, V k * G Y k * and V k = ^i/3|>fe c p'i ) % where Y k is the complementary invariant 
subspace of ker (B* + | J). Furthermore, in the same way as we did for the case k = 0, 
we define the P^ k and P± k , for every k > 1, to be the projections of ker (B* + 1 1) and 
Y k respectively. We also denote V k by 

(5.19) V k *:=n$* 1>k + o(n). 
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Subsequently, substituting (15. 1 8[) into (15.121) and passing to the limit as n \, + , we obtain 
the following equation: 

(5-20) (B* + f = -M k (J2 M=k c P rp), 

under the natural "normalizing" constraint 

(5.21) E c/j = 1. 

|/J|=fc 

Therefore, applying the Fredholm alternative [16] , a unique V fc * G Y" fc * exists if and only if 
the right-hand side of (I5.20p is orthogonal to ker (B* + | /) . Multiplying the right-hand 
side of (I5.20p by if)p, for every \(3\, in the topology of the dual space L 2 , we obtain an 
algebraic system of Mjt + 1 equations and the same number of unknowns, {cp, \(3\ = k} 
and yUi jfc : 

(5.22) <A4(E| /3 |= fc c /3 ^),^> = for all = k, 

which is indeed the Lyapunov-Schmidt branching equation [38]. Through that algebraic 
system we shall ascertain the coefficients of the expansions (I5.10P and, hence, eventually 
the directions of branching, as well as the number of branches. However, a full solution 
of the non-variational algebraic system (I5.22p is a very difficult issue, though we claim 
that the number of branches is expected to be related to the dimension of the eigenspace 
ker(B* + |l). 

In order to obtain the number of possible branches and with the objective of avoiding 
excessive notation, we analyze two typical cases. 

Computations for branching of dipole solutions in 2D . Firstly, we ascertain 
some expressions for those coefficients in the case when = 1, N — 2 and Ml = 2, so 
that, in our notations, {^p}\p\=\ = {^i,^} sucn that ip* = c\ipl + c 2 ?/> 2 . Consequently, 
in this case, we obtain the following algebraic system: 

cSi, h x ) + 2f $ 1: y • V$r) - cxlix,x + c 2 (^i, h 2 ) + £fi (^,y V^ 2 *) = 0, 

(5.23) { Cl (^ 2 , hi) + ^ $ 2 ,y Wj) + c 2 (^2, h 2 ) + sfL $ 2)V . v^*) - c 2/iljl = 0, 
Ci + c 2 = 1, 



where 



hi := -V • [ln(c!-0* + c 2 ^)VA^], h 2 := -V • pn(ci$ + c 2 ^)VA^. 



2J- 



and, Ci, C2 and /xi^ are the coefficients that we want to calculate. Also, ai = Ai is regarded 
as the value of the parameter a denoted by (15. 2p such that 2 are the corresponding 
associated adjoint eigenfunctions. Now, from the third equation we have c 2 = 1 — ci, so 
that substituting it into the first two equations of (14. 5p gives 

(5.24) 



Ni(ci, pn A ) +cif [(iPi,y ■ V^*) - (Vi, y • V^) 



37 



where 



Ni(ci, fj,^) := c^i.hx) + (-01, h 2 ) + ^ (-01 , 2/ • V-0|) - ci(i>i,h 2 ) - ei//i,i, 

-^2(01,^1,1) := Ci(?/>2, frl) + ("02, ^2} + If (-02,2/ • V^ 2 ) - Cl (-02,^2) + Cl^l,! 

represent the nonlinear parts of the algebraic system, with ho and hi just depending on 
Ci in this case. 

To detect solutions for the system ( I5.23|) we apply the Brouwer fixed point theorem to 
( 15.241) (see pQ for further details). Then, we suppose that the values c\ and /ii 5 i are the 
unknowns in a sufficiently big disc Dr(c\, /ti,i), centered in a possible nondegenerate zero 
(cj, /ti i). Therefore, if one of the next two conditions are satisfied 

(^i,2/-V^}-<^i )2 /-W 2 *>^0 or (fo,y-vfe)-{ik,v-vfc)¥:0, 

the nonlinear algebraic system (I5.24p has at least one non-degenerate solution. Note 
that multiplicity results are extremely difficult to obtain. So, to ascertain the number of 
solutions for those nonlinear finite-dimensional algebraic problems is rather complicated. 
However, we expect, and in fact compute it in some cases, that this is somehow related 
to the dimension of the corresponding eigenspace ker (B* + f ij, k > 1. 

Firstly, we calculate the number of solutions for the nonlinear algebraic system (15.231) . 
Integrating by parts the terms in which h\ and h 2 are involved in the first two equations 
and rearranging terms, we arrive at (as usual, all integrals are over ~R N ) 

J Wi • ln(ci^ + c 2 ^)VA(ci^ + c 2 ^)+cif / fay ■ Vfa - ci// M 

+c 2 f fj> iy -Vfe = 0, 



f W> 2 • ln(ci^ + c 2 4*)VA(c 2 ^ + c 2 ^)+cif / i> 2 y ■ - c 2 fii }1 

+c 2 f Ji> 2 yVi>* 2 = 0. 

Then, substituting the third equation with the expression c\ = 1 — c 2 and, hence, putting 
ci^>* + C2-02 — "0* + ("02 — "0i) c 2 into those two equations obtained above yields 

/ V^i • ln(^*+(^ - 4>$)c 2 )VA{4>l + (4* - i){)c 2 ) - At M + C2//1,! 

+f J fay ■ Vfa + c 2 f J fay ■ (V4* - V^) = 0, 



(5.25) 



/ Vfa • ln(^+(4* - ^i> 2 )VA(^ + (fa 2 - fa)c 2 ) - c 2f i 1A 

+f / Al/ • W x * + c 2 f J fay ■ (Vfa 2 - Vfa) = 0. 
Subsequently, adding both equations, we have that 

= /(Wi + V4) • ln(^i + (4* - A*)c 2 ) VA(^* + (fa 2 - fa)c 2 ) 
+f f{fa + • + c 2 ^ /(^i + ^ 2 )y • (V^ - 



:!cS 



Substituting it into the second equation of (15.251) . we find the following equation with the 
single unknown c 2 : 

~4f f$i + fa)y ■ (V^ - V&) + c 2 f(j ^ 2 y ■ Vr 2 - S$ x + 2f 2 )y ■ V&) 

(5.26) +f / ^ 2 y ■ Vfa + j V^ 2 • ln(^ + (4* - ^)c 2 )VA(^ + (4* - ^)c 2 ) 
-c 2 + V^ 2 ) • ln(^ + (4* - ^)c 2 )VA(^ + (4* - ^*)c 2 ) = 0, 

which can be written in the following way: 

(5.27) c\A + c 2 B + C + co(c 2 ) = 5(c 2 ) + u(c 2 ) = 0. 

Here u(c 2 ) can be considered as a perturbation of the quadratic form 5(c 2 ) with the 
coefficients of such a quadratic form defined by 

A _a ffa + fay . - V^), 

B:=f(ffoy. VVr 2 - + 2^ 2 )y ■ V^), 
C:=ffj> 2 y.ViJ>l 

u(c 2 ) := / VV> 2 • lnfrfc + (r 2 - ^)c 2 )VA(^ + (r 2 - ^*)c 2 ) 

-c 2 / (V^i + W> 2 ) ■ Lq(^I + (r 2 - ^)c 2 )VA(^ + (4* - ^*)c 2 ). 

Hence, due to the normalizing constraint ( 15. 21ft . c 2 G [0, 1], solving the quadratic equation 
S"(c 2 ) = yields: 

(i) c 2 = 5(0) = C; 

(ii) c 2 = 1 ^(l) = A + S + C; and 

(iii) differentiating 5 with respect to c 2 , we obtain that $'{c 2 ) = 2c 2 A + B. Then, the 
critical point of the function 5 is c 2 = — £ and its image is 5(c 2 ) = — vr + C. 

Therefore, since we know about the existence of at least one solution, different from 
zero, in this particular case we impose some conditions in order to have at most two 
solutions: 

(a) C{A + B + C) > 0; 

(b) C(-^ + C) <0; and 

(c) < < 1. 

Note that, for — ^ + C = 0, we have just a single solution. 

Consequently, going back again to the equation (15.271) . we need to control somehow 
the perturbation of the quadratic form to maintain the number of solutions. Therefore, 
controlling the possible oscillations of the perturbation u(c 2 ) in such a way that 

||w(c 2 )|| Loo < 5(4), 

we can assure that the number of solutions for (I5.23|) is exactly two (or at most two). 
This is actually the dimension of the kernel for the operator B + ~ I as we conjectured. 
Note that, in general and for large k 3> 1, to solve such multiplicity problems for those 
types of non- variational equations is a rather difficult open problem. 
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Branching computations for 



2. Subsequently, we shall extend those results 



for the case in which the dimension of the eigenspace is greater than 1. Again the calculus 
are rather tedious. For that reason we find it easier to make such calculations for the 
particular case when \/3\ = 2 and = 3 (N = 2), so that {ij)p}\p\=2 — {^ij i>2-> V^s) 
stands for a basis of the eigenspace ker (B* + | J), with k = 2 and A& = — | as the 
associated eigenvalue. Observe that afc(O) = A^. 

Thus, in this case, performing in a similar way as was done for (I5.23p with ip*, = 
Clip* + C21P2 + c&l'ti we arrive at the following algebraic system: 



(5.28) 



Cl(V>l,/ll) +C 2 (^l,/i 2 ) 



+^(4,y-W 3 *>-c 1 // 1 , 2 = 0, 
ci A, M + c 2 (4, ^> + c 2 (^ 2 , /i 3 ) + (^2, y ■ Wj) + (^ 2 , y • V</> 2 *> 

ci(4, fei> + c 2 (^ 3 , fo) + c 2 (^ 3 , ^3) + ^ (4, y ■ Wj) + asa y ■ V^) 



(^3,J/ • V$J) - C 3 /Xl,2 = 0, 



Ci + c 2 + c 3 



where 



/ii := -V • [ln(c!<0* + c 2 tp* 2 + c 3 4*)VA^], /i 2 := -V • pn(ci$ + c 2 ^ + c 3 ^)VA^ 



2J: 



and /iq 



-V • (ln( Cl ^ + c 2 ^ + c 3 ^)VA^ 



3J : 



and ci, c 2 , c 3 , and /ii j2 are unknowns. Here, , 0i,'0 2 ,V> 3 represent the eigenfunctions as- 
sociated with the eigenvalue A 2 = ct 2 (0) and t/>|, ^3 are the corresponding adjoint 
eigenfunctions, which are associated with the same eigenvalue A 2 . 

As for the case \/3\ = 1, the application of the Brouwer fixed point theorem and the 
topological degree provide us with the existence of a non-degenerate solution for the 
nonlinear algebraic system ( 15.28)) under certain conditions. 

Furthermore, in the subsequent analysis, we shall show a possible way to ascertain 
the number of solutions of the nonlinear algebraic system (15.281) . Obviously, since the 
dimension of the eigenspace is bigger than the corresponding one in the case \/3\ = 1, the 
difficulty in obtaining multiplicity results increases. 

We proceed as in the previous case. Firstly, we integrate by parts those terms in which 
hi, h 2 , and h 3 are involved. After rearranging terms, this yields 

/ W>i • ln(ci^ + c 2 </> 2 * + c 3 ^f)V 'A{ci4>* + c 2 ^2 + c 3 ^|)+cif / ipiy ■ V4* - cm, 2 

+c 2 f / fay ■ Vr 2 + c 3 f / fay ■ W| = 0; 



/ W> 2 • ln(ci^ + c 2 fa 2 + c 3 ^)VA(c 2 ^ + c 2 fa 2 + c 3 ^)+cif / fay • Wi - C 2 fMl, 2 

+c 2 f J fay ■ Vfa 2 + c 3 f / fay ■ Vfa = 0; 
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/ V4 • Hci^l + c 2 r 2 + c 3 4*)VA(c 2 ^ 1 * + + c 3 tP* 3 )+ Cl f J 4> 3 y • v4* - C3//1.2 

+c 2 f / ^2/ • W> 2 * + c 3 f / j> 3 y ■ V4* = 0. 

Next, by the fourth equation in (I5.28p . we have that c\ = 1 — c 2 — c 3 . Then, setting 

CiV>l + C 2 4 + C 3 4 = Ipt + C 2 (?p2 ~ 4) + C 3 (4 ~ 4) 

and substituting this into three equations above yields a nonlinear algebraic system: 

/ v4 ■ H#i+$2 - k)<* + (4 - 4*)c 3 )va(4* + (r 2 - r x ) C2 + (4 - 4) C3 ) 

-A*i,2 + c 2 //i )2 + c 3 /ii j2 + sa J . v4 

+ f My • ((V4 - V^)c 2 + (V4 - V4)c 3 ) = 0; 

/ v4 • HkHk - 4*)c 2 + (4 - ^)c 3 )va(^ + (4* - r x ) C2 + (4 - 4) C3 ) 

-c 2 //i, 2 + f S-fay-vfc 

+f J^y • ((V4* - V4)c 2 + (V4 - V4)c 3 ) = 0; 

/ V4 • hi(4 + (4 - 4)c 2 + (4 - 4)c 3 ) VA(4 + (4 - 4)c 2 + (4 - 4)c 3 ) 
(5-29) -W, 2 + ?/^'V^ 

+f / 42/ • ((V4 - V4)c 2 + (V4 - V4)c 3 ) = 0. 
Now, adding the first equation of (15.291) to the other two ones, we have that 

J(v4 + V4) • ln(4 + (4 - 4)c 2 + (4 - ^)c 3 )VA(^ + (4* - 4)c 2 + (4* - 4)c 3 ) 

-// lj2 + c 3 // 1>2 + f f(4 + 4)2/ • Vfc 
+f J(4 + • ((V4 - V4)c 2 + (V4 - V4)c 3 ) = 0, 

J(v4 + V4) • ln(4 + (4 - 4)^2 + (4 - 4)c 3 )VA(4 + (4 - fc)c* + (4 - 4)c 3 ) 

-/ii, 2 + c 2 /ii )2 + f J(4 + 4)2/ • v4 

+f /(4 + 4)2/ • ((v4 - v4)c 2 + (v4 - v4) C3 ) = o. 

Subsequently, subtracting those equations yields 

/ii, 2 = ^ [ /( v4 - v4) • In vl/* VAvT + f /(4 - 4) y ■ v4 

+f J(4 - 4)y • ((v4 - v4)c 2 + (v4 - v4*)c 3 )] , 
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where ** = ^* + — 4) c 2 + (V> 3 — -0*)c 3 . Thus, substituting it into f)5.29p . we arrive 
at the following system, with c 2 and c 3 as the unknowns: 

c 3 / (Wi - W2+ Ws) • In * VA^* -c 2 J (Wi + W2 - Ws) • In #VA^* 
+ / ( W2 - Ws) • In + f /(& - A)?/ • Wi 
+c 2 J [/(A - 4)y • V(4* - 2&) - /&y • Wi] 
+c 3 f [/(A - fa)y ■ V(4 - 2-0*) -jfay Wi*] 
+e 2 c 3 f [/^-(Wf-Wa*) 
+ /(^ - fa)y ■ (2V^ - W 2 * - W 3 *))] 

+cif /(A - 4 + 4)y ■ (W 3 - W?) 



—c 



2a 2 
2 4 



+ 4 - 4)2/ • (W 2 * - W?) = 0, 



c 3 / V4 • In VA^* - c 2 / Ws • In ^* VA#* + c 3 ^f / fay • Wi 

-c 2 f / fay ■ Vfa + c 3 f / 4?/ • (( W 2 * - Wi> 2 + ( W 3 - Wi)c 3 ) 

-c 2 f /^y • ((W 2 * - Wi> 2 + (W 3 * - Wi> 3 ) = 0. 

These can be re-written in the following form: 

. ^l(c2,C 3 ) +Wi(c 2 ,C 3 ) =A 1 c\ + Bxcl + 0x02 + 0^ + E1C2C3 + wi(c 2 , c 3 ) = 0, 

(5.30) 

3i(c 2 , c 3 ) + w 2 (c 2 , c 3 ) = A 2 c 2 + 5 2 c 3 + C 2 c 2 + D 2 c 3 + E 2 c 2 c 3 + cu 2 (c 2 , c 3 ) = 0, 



where w 1 (c 2 , c 3 ):= c 3 / (Wi - W2 + Ws) • In tf*VAtf* 
-c 2 / (Wi + W 2 - W 3 ) • In ^VAtf* 
+ /(W2 - WO • In** - f /(4 - 4)2/ • Wi, 

and u; 2 (c 2 ,c 3 ):= c 3 / W2 • ln^*VA^* - c 2 / W 3 • ln^*VA#* 
are the perturbations of the quadratic polynomials 

&(c 2 , c 3 ) := AjC2 + SjCg + CiC 2 + Ac 3 + £;C 2 c 3 , with 2 = 1,2. 
The system (I5.30P can be re-written in a matrix form with two quadratic forms involved: 

fa c 3 )Pi(;) + Qi(3+Fi = o, 

(c 2 c 3 )P 1 g)+Q 1 g)+F 2 = 0, 

where the matrices Pj and Q., of the quadratic forms with j = 1,2 have the corresponding 
coefficients Aj to Ej as entries, plus the perturbations of the quadratic forms denoted, 
under this notation, by Fj, with j = 1,2. 

Then, owing to the conic classification, we are able to solve (I5.30P (without the nonlinear 
perturbation) and obtain an estimate for the number of solutions of the original nonlinear 
algebraic system f l5.28[) . 
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Hence, according to the conic classification, we will have the following conditions that 
will provide us with conic section of each equation of the system (I5.30p (without the 
nonlinear perturbation): 

(i) If Bj—AAjCj < 0, the equation represents an ellipse, unless the conic is degenerate, 
for example c\ + c§ + a = for some positive constant a. So, if Aj = Bj and Cj = 
the equation represents a circle; 

(ii) If B? — 4:AjCj = 0, the equation represents a parabola; and 

(iii) If Bj —AAjCj > 0, the equation represents a hyperbola. If we also have Aj + Cj = 
the equation represents a hyperbola (a rectangular one). 

Therefore, taking into account the "normalizing" constraint, the zeros of our system 
will depend on the coefficients we have for the system, so on the eigenf unctions that 
generate each eigenspace ker (B* + | l) . 

Observe that the number of intersections between two conies oscillates from one to four. 
Hence, this will be the possible number of branches that are obtained for our problem. 
However, since the dimension of the eigenspace in this particular case is three, it seems 
that, in this case, we have four branches, so two of them should coincide, though this 
claim remains uncertain. 

Moreover, as was done for the previous case when \/3\ = 1, we need to control the 
oscillations of the perturbation functions in order to maintain the number of solutions. 
Consequently, assuming that 

m(c 2 ,c 3 )|| LOO < $i(c* 2 ,c* 3 ), with i = l,2, 

we conclude that the number of solutions must be between one and four. This again gives 
us an idea of the difficulty of more general multiplicity results. 

6. Final comments 

6.1. A first comment: towards evolutionary completeness. According to |25J , evo- 
lutionary completeness of the nonlinear eigenfunction subsets $~(n) simply means that 
those functions describe all possible types of finite time blow-up asymptotics for solutions 
of the TFE-4 fll.ip in a neighbourhood of any point (xo,to). For nonlinear evolution 
equations, such a completeness is a very difficult question. As far as we know, the evolu- 
tion completeness result proved in [25] for the ID porous medium equation on a bounded 
interval remains the only known such result for essentially quasilinear PDEs (i.e., not a 
perturbed semilinear equation). 

Indeed, for the TFE-4 (11.11) . such a completeness problem is difficult beyond any imag- 
ination. In particular, this would include a full analysis of all the asymptotics of the 
non-stationary quasilinear fourth-order degenerate parabolic flow (12.101) containing no 
monotone, coercive, potential, or order-preserving operators. 

However, our homotopy approach somehow implies certain (but not that strong and 
promising) confidence concerning the evolutionary completeness of $ _ (n) for n > 0: 
the point is that, for n = 0, the eigenfunction set of the Hermite polynomials (13.171) is 
indeed complete and closed in any suitable weighted space, where those notions are now 
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understood as in classic theory of bi-orthogonal polynomials and Riesz bases. So we may 
expect that the evolution completeness for small n > can be "inherited" from those 
brilliant spectral properties available for n = (Section [3]). This is the only issue we are 
aware of and can rely on in this analysis. 

Same speculations apply to the evolutionary completeness of global similarity patterns 
for small n > 0, which is now connected with completeness/closure of eigenfunc- 
tions (I3.15P of B for n = 0; see [18J for proofs. 

6.2. A pessimistic comment. Overall, we must admit that, though we have obtained 
some multiplicity results for not-that-multi-dimensional eigenspaces and have shown cer- 
tain extensions of our techniques, any further rigorous justification seems to be too ex- 
cessive. Indeed, any rigorous results will inevitably require to specify or evaluate with 
sufficient accuracy of those numerical values of various projections given by linear func- 
tional as linear combinations of the Hermite polynomials ( 13.1 7p . In view of a complicated 
nature of non-self-adjoint theory for the spectral pair {B, B*}, this is expected to be en- 
tirely illusive. 

On the other hand, it would be very important to trace our n-bifurcation branches 
of nonlinear eigenvalue problems in both global and blow-up setting by using some more 
general and powerful techniques of nonlinear operator theory. However, no one can expect 
this to be a simply task. We suspect that, in view of principally non- variational structure 
of such nonlinear eigenvalue problems, containing no monotone and/or strongly coercive 
operators, any non-local (in n) sharp results on existence/multiplicity of n-branches will 
not be obtained reasonably soon. 

Therefore, overall, we claim that our n-branching approach, which allowed us to explain 
the occurrence of nonlinear branches from linear eigenfunctions at n = 0, though not being 
fully rigorous, is the only currently available way to detect branching phenomena for 
such nonlinear eigenvalue problems embracing similar classes of non-variational and non- 
monotone operators. It is clear how these homotopy-branching methods can be extended 
to more general and more higher-order quasilinear operators of different types, once a 
parameter homotopy to a proper linear spectral problem for a suitable non-self-adjoint 
pair {B,B*} has been well understood and carefully and rigorously studied. However, 
we warn that the latter linear problem often can be a very difficult one itself; one such 
example of a refined scattering theory for 2m-th order linear Schrodinger operators is 
under attack in [27] . 
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Appendix A. Necessary functional setting for branching at n = + 

Here, we are going to present some justification of the our main branching analysis. Namely, 
we need to deal with expansions such as (jl.20p and/or (jl.2ip . Recall that, using this, we are 
not going to, and in fact cannot, justify rigorously the existence of nonlinear eigenfunctions as 
solutions of (|1.3p , , (|1.4p at least for small n > 0, but just the branching at n = 0, under the 
hypothesis that a proper limit 

(A.l) f(y; n) -> f (y) (= ^(y)) as n 0+ 

exists in a necessary metric to be specified. According to our spectral theory of the non- 
self adjoint pair {B,B*}, here /o denote some eigenfunction ipp, and, in the most simple and 
interesting case, we assume that 

(A.2) Mv)=F(y)(=Mv)), 



where F{y) is the rescaled kernel of the fundamental solution (|3.3p of the bi-harmonic operator. 
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Thus, we need to check under which extra assumptions on (|A.ip . the following limit takes 
place, in the weak sense, 

(a.3) \m^±^i n \ foiy) i 

where the right-hand side is assumed to be well defined (bounded) a.e. First of all, it is obvious 
that such a convergence crucially depends on the structure of zeros of the limit functions fo(y), 
which is easy to demonstrate: 

Example: a non-transversal zero. Let fo(y) have a non-transversal zero at, say, y = 0~ 
(the interface point), and 



(A.4) /( y;n )=e^ for y < => = ^ -» oo. 

Actually, this means that 

(A. 5) f(y; n) -> = f {y) as n -> + for all y < 0, 

i.e., y = is not a transversal zero of fo(y). Then the limit (\A.3\i makes no sense and the 



branching analysis at re = does not apply at all. 

Fortunately, such a situation cannot occur for the analytic kernel F{y) and all its derivatives, 
representing other eigenfunctions. Of course, we cannot guarantee that non-transversal zeros of 
F{y) cannot occur at all. They can, but with a lower probability as for any analytic function. 
However, we do know that such non-transversal zeros are always isolated and cannot concentrate 
on a given surface in R . Therefore, on any compact subset such non-transversal zero surfaces 
have zero measure. However, this is not sufficient and an extra rough estimate would be useful. 

Evidently, (|A.3P is violated in the pointwise sense on a bad set of points b*(n) such that 

(A.6) \f(y; n)\ « |/ (y)| <. e"« for all re > small. 

Then, in this worst case, 

(A.7) l/o(y)l"-i „ _ 1 _> oo as n ^ +. 

Assume that a non-transversal (a multiple) zero again occur at y = and the ID behaviour is 
as follows: 

(A.8) fo(y)~y\ where k = 2,3,4, ... . 

Then, in the weak sense, the integral representation of (|A.3j) will provide us with the "bad" 
( "worst" ) discrepancy of the order 

(A. 9) ~ i meas&*(n) ~ ^ e - ^ — > as re — )• + 

for any finite multiplicity of the zero at y = 0. It is clear that any use of the M^-geometry of such 
multiple zeros cannot affect the non-analytic exponential term in (|A.9h and the convergence. 
We complete our discussion as follows: 

Proposition A.l. Let ()A.1|) hold uniformly on compact subset, where the limit function fo(y) 
satisfy the above assumption of a. a. transversal zeros. Then (|A.3P holds in the weak sense. 
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Finally, let us also formally note that, in ()1.24|) on the bad set 6*(n), we have the following: 

(A.10) ((VA)-^^ VA/)| ~ ((VAr^IVA/)! ~ \lf(y)\ « ± e"£ -> 

as n — > + . This confirms that the convergence (|1.24p takes place a.e., provided that the zero set 
of fo(y) has zero measure only on any compact subset in R , i.e., the analyticity is not required 
(Sard's theorem for C p functions in R^ with any p > 1 may be used). Of course, this is just a 
rough estimate and further study is needed. 
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